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SUMMARY 


A series  of  experiments  and  analyses  were  carried  out  on  the  cyclic 
crack  growth  behavior  of  center- cracked  cruciform  specimens  under  various 
biaxial  loading  conditions  (-1.5  < ox/o  ^ 1.75  for  constant  amplitude. tests, 

0.5  < ax/Oy  <0.5  for  periodic  single  overload  and  spectrum  load  tests). 

The  results  may  be  summarized  as  follows: 


1.  For  cracks  perpendicular  to  a the  crack  grew  straight  except 

for  a < o • 
x — y 


2.  Elastic  K factors  were  obtained  for  both  straight  and  curved  cracks 
and  were  adequate  for  correlating  the  biaxial  crack  growth  rate 
data. 


3.  The  direction  of  crack  growth,  the  crack  tip  stress  intensity  and 
the  crack  growth  rate  were  controlled  by  the  most  dominant  stress 
component  in  a multiaxial  stress  field.  For  a straight  crack,  the 
effect  of  ox  on  constant  amplitude  crack  growth  rate  was  negligible 
provided  the  crack  tip  stress  intensity  factors,  in  a given  biaxial 
stress  field,  were  properly  computed.  Crack  growth  rates  for  the 
curved  and  the  straight  cracks  were  the  same  at  a given  stress 
intensity  level. 

4.  In  out- of- phase  loading  conditions,  the  crack  growth  rate  and  crack 
growth  directions  (for  ox  < a ) were  the  same  as  those  in  the 
in-phase  loading  condition. 

5.  Crack  tip  plastic  zone  size  variations  with  biaxial  ratio  appeared 
to  have  no  effect  on  constant  amplitude  crack  growth  rate. 

6.  At  lower  stress  intensity  levels,  the  analytically  determined  crack 
tip  plastic  zone  sizes  qualitatively  correlated  with  those  directly 
measured  from  the  test  coupons  (by  using  an  interferometry  micro- 
scope). At  higher  K levels,  the  validity  of  the  interferometry 
measured  plastic  zone  sizes  was  questionable. 
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7.  As  for  the  variable  amplitude  tests,  it  was  very  evident  that  cracks 
grew  faster  at  positive  biaxial  stress  states  but  slower  at  negative 
biaxial  loading  conditions  compared  to  uniaxial  loading  conditions. 
The  crack  growth  retardation  behavior  were  similar  in  both  uniaxial 
and  biaxial  states  of  stresses;  i.e.,  the  delay  cycles  were 
controlled  by  the  load  interactions.  Therefore,  existing  retardation 
models  could  be  used  to  predict  delay  cycles  in  both  uniaxial  and 
biaxial  stress  fields  provided  that  the  crack  tip  plastic  zone  sizes 
in  a given  stress  state,  were  accurately  determined. 
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SECTION  I 


INTRODUCTION 


Fracture  mechanics  techniques  are  currently  being  used  to  perform  damage 
tolerant  analysis  on  many  types  of  aircraft  structural  components.  Because 
the  problems  of  multiaxial  loading  are  very  common  in  aircraft  airframe  and 
engine  components,  it  is  important  to  evaluate  and  quantify  multiaxial  effects 
in  order  to  improve  the  crack  propagation  prediction  capability  for  design 
purposes. 

Consider  that  a plate,  containing  a through- the- thickness  crack,  is 
subjected  to  a biaxial  stress  field.  One  of  the  stress  components  is  acting 
perpendicular  to  the  crack  and  another  component  parallel  to  the  crack.  For 
problems  that  are  solved  using  purely  elastic  formulation  [l,  2],  the  crack 
tip  stress  intensity,  K^,  in  the  opening  mode,  is  theoretically  not  affected 
by  the  lateral  stress  component.  On  the  other  hand,  it  can  be  shown  by 
elastic-plastic  analysis  [2-8]  that  the  size  of  a crack  tip  plastic  zone 
varies  with  biaxial  loading  conditions.  Therefore,  one  may  speculate  that 
the  crack  tip  stress  intensity  (K)  as  well  as  the  crack  growth  rate  (da/dN) 
will  also  be  influenced  by  the  presence  of  the  lateral  stress  component. 
Experimental  data  concerning  the  biaxial  loading  effect  on  cyclic  crack 
growth  rate,  residual  strength,  and  the  direction  of  crack  growth,  are 
available  [9-19].  However,  the  results  of  these  few  investigations  have  been 
inconsistent,  and  none  of  them  contained  enough  data  to  offer  conclusive 
evidence  to  support,  or  to  correlate  with,  the  existing  theories. 

The  objective  of  this  investigation  was  to  systematically  evaluate  the 
effect  of  the  biaxial  stress  field  on  cyclic  crack  growth  rate  behavior.  An 
analytical/experimental  study  program  was  designed  to  answer  the  following 
questions: 

1.  Do  biaxial  stresses  affect  crack  tip  stresses  intensity,  cyclic  crack 
growth  rate,  crack  growth  direction,  or  crack  tip  plastic  zone  size? 

2.  Is  an  elastic  - K adequate  for  correlating  the  biaxial  da/dN  data? 

3.  Is  the  crack  tip  plastic  zone  important  in  the  mechanics  of  crack 
growth? 
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SECTION  II 


STATEMENT  OF  WORK 


An  exploratory  research  program  was  conducted  to  systematically 
evaluate  the  effects  of  biaxial  stress  ratio  on  constant  amplitude  and 
variable  amplitude  fatigue  crack  growth  rates.  The  scope  of  the  program 
included  the  following: 

1.  The  specimen  design  was  optimized  by  finite  element  analysis.  The 
analytic  results  also  determined  the  load  and  stress  relationship  on 
the  optimized  specimen  geometry. 

2.  A strain  survey  was  conducted  to  verify  the  analytic  load-stress 
relationship. 

3.  A crack-finite-element  technique,  and  other  appropriate  means,  were 
used  to  determine  elastic  crack  tip  stress  intensity  factors  for  a 
crack  in  the  selected  specimen  geometry.  Crack  geometries  included 
the  straight  crack,  inclined  crack,  crack(s)  at  a hole,  and  curved 
crack. 

4.  A total  of  118  specimens  made  of  7075-T7351  and  2024-T351  aluminum 
alloys  were  tested  under  various  biaxial  loading  conditions. 

A vast  amount  of  cyclic  crack  growth  rate  and  fracture  data  was 
developed.  The  scope  of  the  test  program  is  outlined  in  the 
following: 

A.  Evaluation  of  the  basic  biaxial  effects  by  conducting  constant 
amplitude  fatigue  crack  growth  tests  at  various  biaxial  stress 
ratios  (-1.5<  B < 1.75),  at  R =0.1. 
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B.  Evaluation  of  the  basic  biaxial  effects  by  conducting  constant 
amplitude  fatigue  crack  growth  tests  at  R = 0.7  (i.e.,  evaluation 
of  the  R effect). 

C.  Evaluation  of  the  biaxial  ratio  effects  on  fatigue  crack(s) 
coming  out  from  a hole  configuration. 

D.  Evaluation  of  crack  growth  retardation  behavior  under  biaxial, 
variable  loading  conditions  (with  -0.5  < B < +0.5). 

E.  Evaluation  of  the  biaxial  stress  effect  on  constant  amplitude 
fatigue  crack  growth  rate  under  out-of-phase  loading  conditions. 

F.  Miscellaneous  tests  including  sustained  load  tests,  angle  crack 
in  a biaxial  stress  field,  and  fracture  tests. 

5.  The  "loaded”  crack  tip  plastic  zone  sizes,  at  various  biaxial  stress 
ratios,  were  determined  by  using  the  elastic-plastic  finite 
element  modeling  technique.  The  "loaded"  and  the  "unloaded"  crack 
tip  plastic  zone  sizes  were  also  measured  from  the  surface  of  a 
specimen  (during  a cyclic  loading  test)  by  using  an  interferometry 
microscope . 

6.  The  experimental  data  were  analyzed  to  determine  whether  or  not 
biaxial  stress  would  affect  the  fatigue  crack  growth  behavior.  A 
computer  program  was  prepared  to  compute  the  crack  growth  rate, 
as  well  as  the  crack  growth  direction,  under  constant  amplitude 
biaxial  loading.  Variable  amplitude  fatigue  crack  growth  rate  data 
were  correlated  with  the  modified  Willenborg  model  and  a new  crack 
growth  retardation  model  was  developed. 

7.  Based  on  the  results  of  this  research  program,  guidelines  and 
procedures  for  predicting  safe-crack-growth- life  in  a biaxial  stress 
field  were  derived. 
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SECTION  III 


BIAXIAL  TESTING  MACHINE 

A biaxial  test  machine,  shown  in  Figures  1 to  3,  was  made  specifically 
for  testing  of  the  specimens  described  in  Section  4.  The  unit  consists  of 
three  hydraulic  actuators  mounted  on  a horizontal  frame  made  of  aluminum 
channel  (provisions  for  a fourth  hydraulic  actuator  were  integrated  into 
the  test  machine  design).  Two  of  the  actuators  are  plus  and  minus  55  kip 
capacity  mounted  opposing  each  other  (this  direction  is  designated  the  "X" 
direction).  The  third  actuator  is  plus  and  minus  77  kip  capacity  and  is 
mounted  perpendicular  to  the  other  two  actuators  in  the  "Y"  direction.  No 
actuator  (presently)  opposes  the  plus  and  minus  77  kip  actuator;  instead, 
this  end  is  fixed  (similar  to  a typical  crosslead).  Each  actuator  mount 
is  attached  to  the  test  frame  through  pin-bearing  arrangements  (upper  and 
lower)  which  allow  the  actuators  to  swivel  (in  a horizontal  plane).  This 
was  required  to  prevent  possible  actuator  damage  from  side  loading  (which 
could  occur  during  specimen  failure,  for  example).  The  bearing-pin  fit  is 
approximately  a 0.001  inch  clearance,  which  minimizes  load  train  stop. 

Each  actuator  is  terminated  into  a large  device  (spiral  locking  rings  are 
used  for  proper  prestressing  of  threaded  adapters).  The  device  holds  a 
second  pin-bearing  arrangement  (similar  to  the  actuator  pin-bearing  arrange- 
ment); this  pin  serves  as  the  loading  pin,  which  the  grips  contact.  The 
specimen  loading  arms  are  "sandwiched"  between  these  grips  in  a seven-bolt 
hole  pattern.  The  top  cover  plate  is  removable  for  specimen  set-up  and  is 
locked  down  after  such.  The  fixed  end  of  the  test  frame  is  constrained  so 
that  the  device  cannot  swivel.  This  is  required  for  tests  involving  com- 
pression loading. 

Three  load  cells  are  used  in  the  test  frame.  These  are  mounted  to  the 
device  on  both  55  kip  ends  and  on  the  fixed  end,  (i.e.,  each  load  cell  is 
mounted  on  the  device  opposing  the  corresponding  actuators).  All  load  cells 
are  plus  and  minus  100  kip  fatigue  rated  GSE  load  cells.  Each  actuator  is 
also  equipped  with  an  LVDT,  and  is  controlled  by  a 15  gpm  Moog  servo-valve. 
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55  KIP 


Biaxial  Loading  Machine 


Close  Up  Of  The  Biaxial  Loading  Machine 


Figure  3.  Control  Unit 


The  three  actuators  are  controlled  from  a console  consisting  of  three 
MTS  Model  406  servo  controllers  and  one  MTS  Model  436  control  unit.  Dual 
MTS  Model  430  digital  indicators  were  used  for  monitoring  of  load  and 
stroke  signals.  An  MTS  Model  410  digital  function  generator  was  used  for 
commanding  of  constant  amplitude  and  static  (strain  survey)  tests.  An 
automated  control  unit  was  used  for  commanding  of  spectrum  tests. 

For  all  tests,  the  mode  of  control  of  the  system  is  described  as  follows. 
The  77  kip  and  one  of  the  55  kip  actuators  were  run  in  load  control 
commanded  by  the  servo-controller  (Model  406)  load  control.  The  second 
55  kip  actuator  was  run  in  stroke  control.  The  command  signal  to  this 
actuator  was  the  stroke  feedback  from  the  opposing  55  kip  actuator.  This 
method  of  control  was  employed  so  that  the  two  55  kip  actuators  displaced 
in  continuous  phase;  therefore,  the  specimen  centerpoint  (X-direction)  was 
held  fixed  and  thus  maintained  proper  specimen  positioning  for  loading  in 
the  Y-direction.  Note,  however,  that  the  specimen  centerpoint  in  the 
Y-direction  is  not  so  constrained  (since  there  is  no  fourth  actuator). 
Therefore,  some  off-center  line  loading  was  imposed  in  the  X-direction. 

From  check  outs,  it  was  determined  that  the  lateral  displacement  of  the 
55  kip  actuator  was  approximately  0.001  inch  per  kip  (77  kip  load)  at  the 
specimen  centerpoint.  Considering  the  distance  from  specimen  center  to 
actuator  swivel  point  (about  27  inches),  the  angular  misalignment  is  quite 
small  and  is  considered  insignificant. 

A typical  example  of  test  set-up  is  described  below.  The  precracked 
test  specimen  was  mounted  in  the  biaxial  test  frame  using  low  pressure  load 
control.  Small  tensile  loads  (approximately  500  pounds)  were  applied  in 
both  directions.  To  set  the  specimen  position  in  the  X-direction,  the 
55  kip  actuator  in  stroke  control  was  positioned  so  that  the  77  kip  actuator 
had  minimum  stroke  displacement.  This  minimum  occurred  when  the  77  kip 
actuator  was  in  line  with  the  opposing  fixture  device.  Final  bolt  tightening 
was  performed  and  the  loads  were  returned  to  zero  level. 
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The  load  cycling  was  intitiated  and  the  spans  (amplitude  control 
potentiometers)  were  opened  in  small  increments  (alternating  between  X and 
Y actuators).  The  two  digital  peak  indicators  monitoring  load  were  used 
to  indicate  actual  load  response.  Five  adjustments  on  the  spans  were 
made  until  the  desired  load  peaks  were  accurate  to  one  percent. 

Phase  angle  is  optimized  by  adjusting  servo-controller  gains  and  compar- 
ing X and  Y load/ time  traces  on  a storage  oscilloscope.  Such  gains  are 
"tweeked"  in  until  no  discernible  phase  lag/lead  is  noticed  (this  applies  to 
0°  phase  testing;  for  180°  phase  the  feedback  signals  for  one  of  the  load 
cells  is  inverted  to  permit  a comparison).  Finally,  the  spans  usually 
require  a slight  readjustment  after  gain  adjustment  to  maintain  load  peaks 
accurate  to  better  than  plus  and  minus  one  percent. 

This  complete  procedure  is  accomplished  in  a matter  of  a few  minutes  at 
most  (as  more  experience  was  gained  this  time  interval  was  reduced 
considerably) . 
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SECTION  IV 


SPECIMEN  CONFIGURATION 


There  are  many  types  of  specimens  that  we  can  use  to  accomplish  a biaxial 
loading  condition.  For  example,  Pook  and  Holmes  M used  a flat  cruciform 
specimen  containing  longitudinal  slots  in  the  loading  arms;  Beck  M used, 
a very  large  square  sheet  and  loaded  the  sheet  through  many  little  straps 
attached  around  the  sheet  edges.  The  criteria  for  designing  a specimen 
configuration  to  fulfill  all  the  objectives  in  the  present  study  were: 

1.  The  specimen  should  be  capable  of  taking  compression  load. 

2.  The  specimen  should  be  designed  to  avoid  fatigue  damage  at  the  grip 
or  in  any  area  other  than  that  containing  the  crack. 

3.  The  size  of  the  specimen  should  be  large  enough  to  minimize  boundary 
effects  on  crack  tip  stress  intensity;  but  it  should  not  be  too  large, 
so  that  the  required  load  levels  can  be  kept  within  the  capacity  of 
the  testing  machine. 

4.  The  stress  distribution  across  the  specimen  width  should  be  fairly 
unif  orm. 

5.  The  specimen  configuration  should  be  simple  in  order  to  minimize 
machining  costs. 

A cruciform  specimen  configuration  was  selected  for  generating  biaxially 
loaded  crack  growth  rate  data.  Generally,  the  specimen  had  an  overall 
length  of  597  mm  (23.5  inches)  including  grip  areas  at  each  end  of  the  loading 
arms.  It  also  had  a thinner  region,  152.4  mm  (6  inches)  in  diameter,  in  the 
center  of  the  specimen.  An  overall  view  of  the  specimen  is  shown  in  Figure  4. 
Figure  5 is  a close  up  photograph  of  the  center  section.  Loading  conditions 
and  dimensions  of  the  specimen  are  shown  in  Figure  6. 
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Figure  6.  Finite  Element  Model  For  One-Quarter  of  a Cruciform  Specimen 


It  was  determined  that  the  thickness  of  the  center  region  (t^)  and  the 
thickness  and  width  of  the  loading  arms  ( t 2 and  W)  were  the  three  primary 
variables  affecting  the  stress  distribution.  A 17.7  mm  (0.5  inch)  thick 
loading  arm  was  selected  for  t2  to  eliminate  one  of  the  three  variables  and 
also  to  minimize  material  and  machining  costs.  The  t^  and  W dimensions  were 
optimized  by  conducting  stress  analysis  on  a dummy  panel  configuration 
(without  crack).  Stress  distributions  across  the  thin  section  were  determined 
by  using  the  NASTRAN  computer  program.  Figure  6 shows  the  finite  element 
model  representing  one  quarter  of  the  cruciform  specimen.  Here  t^  = 4.57  mm 
(0.18  inch),  W = 17.78  cm  (7  inches)  and  t2  = 17.7  mm  (0.5  inch).  Also 
shown  in  Figure  6,  are  three  rings  of  triangular  elements  of  different 
intermediate  thicknesses  (t^»  t^  and  t^)  to  simulate  the  curvature  connecting 
t^  and  t2* 


The  analytical  results  are  presented  by  the  curves  of  Figures  7 and  8 . 
In  Figure  7,  the  load  and  stress  relationship  at  the  center  of  the  specimen 


plotted  as  functions  of  P^  to  P^  ratio.  The  load  and  stress  relationship 
(the  NASTRAN  lines  in  Figure  7)  can  be  represented  by  the  following  equations. 


cr 

y 


a 

x 


(1) 

(2) 


For  an  actual  test,  the  required  P^  and  P values  corresponding  to  any 
desirable  ox  and  a combinations  can  be  determined  by  solving  Equations  1 
and  2. 


In  Figure  8,  stress  distributions  along  the  X-axis  of  the  cruciform 
specimen  are  presented.  Since  the  specimen  was  symmetrical  about  its  center 
lines,  the  magnitudes  of  crx  and  a are  plotted  as  functions  of  a.  The  stress 
distributions  corresponding  to  many  Px/Py  ratios  were  determined;  however, 
only  four  typical  examples  are  shown  here.  A complete  set  of  stress 
distribution  curves  is  presented  in  Volume  II  of  this  report. 


14 


(NASTRAN)  \<t,  (NASTRAN) 


Figure  8.  Stress  Gradient  along  the  X-Axis  of  a Cruciform  Specimen 


SECTION  V 


STRESS  SURVEY 


One  cruciform  specimen  of  each  material  was  instrumented  with  rectangular 
rosettes  along  the  X-axis  of  the  specimen.  The  specimens  (without  crack  or 
cut,  or  elox)  were  loaded  in  the  horizontal  biaxial  loading  machine  under 
various  biaxial  load  ratios  (e.g.,  0,  + 0.5,  + 1.0,  + 2.0,  etc.).  Stress 
distributions  along  the  X-axis  of  the  specimen  were  determined  by  strain 
measurements  at  five  locations,  one  pair  of  back-to-back  rosettes  at  the 
center  of  the  specimen  and  two  rosettes  on  each  side  of  the  center  covering 
a 69.85  mm  (2.75  inch)  radius.  Two  of  the  side  locations  also  had  back-to- 
back  rosettes  (see  Figure  9).  A guide  was  used  to  prevent  out-of-phase 
plane  buckling  under  compression  loads. 

At  each  loading  condition,  up  to  four  load  levels  were  applied  and  two 
readings  were  taken  at  each  load  level.  The  specimen  was  placed  in  the 
machine  at  two  orientations.  One  set  of  the  strain  gauge  data  was  taken 
while  the  Y-axis  of  the  specimen  was  lined  up  with  the  244,750  N (55  kip) 
load  cells  and  another  set  of  strain  gauge  data  was  taken  while  the  Y-axis 
of  the  specimen  was  lined  up  with  the  342,650  N (77  kip)  load  cells  (i.e., 
the  specimen  was  rotated  90°).  Typical  experimental  data  are  plotted  on 
Figures  7 and  8.  It  is  observed  in  these  figures  that  the  correlations 
between  strain  gauge  results  and  the  NASTRAN  finite  element  analysis  results 
are  very  good.  Note  that,  in  these  figures,  each  data  point  represents  an 
average  of  two  readings.  Occasionally,  there  is  a number  adjacent  to  a data 
point  indicating  that  more  than  one  data  point  was  superimposed  on  another; 
e.g.,  the  number  four  implies  that  the  data  point  represents  an  average  of 
eight  measurements  (four  load  levels  and  two  reading  per  each  load  level). 
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The  solid  and  the  open  symbols  in  Figure  8 indicate  the  results  from  back- 
to-back  gauges.  Since  the  results  for  the  back-to-back  gauges  were  almost 
identical,  only  one  side  of  the  results  is  presented  in  Figure  7.  It  is 
significant  to  note  that  experimental  data  (a  complete  set  of  the  data  is 
presented  in  Volume  II  of  this  report)  showed  that  the  load  response 
characteristics  of  the  cruciform  specimen  were  not  affected  by  the  position 
of  the  specimen;  i.e.,  whether  the  Y-axis  of  the  specimen  was  placed  in  line 
with  the  342,650  N(77  kip)  or  the  244,750  N(55  kip)  load  cells,  the  strain 
gauge  results  were  identical. 
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SECTION  VI 


STRESS  INTENSITY  FACTORS 


6.1  Straight  Crack 

In  the  theory  of  linear  elastic  fracture  mechanics,  crack  tip  stress 
intensity  can  be  expressed  as 


K = CTyjTTa  * F (a) 


where  ay  is  the  gross  area  stress  normal  to  the  crack.  In  case  of  a center 
cracked  specimen  (CCT) , O’  will  be  the  far  field  uniform  stress  and  F(a) 
accounts  for  the  boundary  conditions.  According  to  [20] , 


F(a)  = [l  - 


(A) 


In  case  of  a cruciform  specimen,  for  a given  ratio  of  biaxial  loads, 
there  is  a pair  of  stress  components,  ax  and  ay » at  every  point  along  a 
predetermined  crack  plane.  In  this  case,  the  crack  plane  will  be  the  X-axis 
in  Figure  6 . As  postulated  in  w.  the  elastic  K value  for  a given  crack 
length  in  a biaxial  state  of  stress  should  be  the  same  as  in  the  uniaxial 
loading  condition.  In  other  words,  the  K-expression  of  Equation  3 is 
applicable  to  the  cruciform  specimen  except  that  O’ y would  be  the  (reference) 
stress  in  the  center  of  the  uncracked  specimen  and  F(a)  would  be  a function  of 
the  boundary  conditions  and  the  stress  gradient  of  CTy  along  the  X-axis. 

A finite  element  analysis  of  the  cruciform  specimen  with  cracks  was 
conducted  using  the  finite  element  model  of  Figure  6 to  determine  elastic 
K values.  A special  "crack  tip"  element,  originally  developed  by  Tong, 
et  al,  [.21]  was  incorporated  into  the  NASTRAN.  In  each  case  analyzed, 
e.g  , each  crack  length,  a special  element  was  placed  in  the  general 
finite  element  model  occupying  a region  representative  of  the  predeter- 
mined crack  tip  location,  and  the  elements  representing  the  crack 
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were  freed  from  the  boundary  restrictions.  Eight  specimens  with  crack 
lengths  (a  = 6.35  to  69.85  mm)  were  loaded  to  various  biaxial  loading  ratios 
with  cr  = 12  ksi.  The  results  are  plotted  in  Figure  10.  Several  K values 
calculated  from  Equation  3 with  F(a)  = 1.0  are  also  plotted  in  Figure  10  for 
comparison.  It  is  seen  that  the  effects  of  loading  conditions  on  elastic  K 
values  are  negligible  and  that  the  cracked  cruciform  specimen  behaves 
actually  like  an  infinite  sheet  especially  at  positive  biaxial  loading 
conditions.  It  is  even  more  important  to  note  that  the  apparent  deviations 
in  K,  for  a crack  length  greater  than  38.1  mm  (1.5  inches),  were  mainly  due 
to  the  effect  of  specimen  geometry  rather  than  the  effect  of  biaxial  loading 
ratios.  The  hypothesis  is  substantiated  by  the  fact  that  the  K values  for 
long  cracks  under  negative  o~x  loads  were  actually  lower  than  those  under 
positive  crx  loads.  Comparing  Figures  8A  and  8B  to  Figures  8C  and  8D,  it  is 
evident  that  the  tension-compression  loading  cases  exhibited  more  reductions 
in  the  cr  stresses  in  the  area  near  the  rim. 

y 


6.2  Angle  Crack 

Table  1 shows  the  cracked  finite  element  analysis  results  for  a 45° 

crack  in  a cruciform  specimen.  Three  loading  cases  were  analyzed;  i.e., 

<T  = 0,  a = a and  a - -a  . In  all  three  cases,  the  input  stress  for  a 
x x y x y y 

was  12  ksi.  Due  to  the  unsymmetric  nature  of  the  problem,  a full  specimen 
model  was  used  (see  Figure  11).  Stress  intensities  at  both  ends  of  the  crack 
were  computed.  As  shown  in  Table  1 , the  calculated  K- values  were  almost 
the  same  at  each  crack  tip.  Theoretically,  in  an  infinite  sheet,  the  <Jx  = CF 
loading  condition  will  be  pure  Mode  1 and  the  <7^  = -q  loading  condition 
will  be  pure  Mode  2.  As  can  be  seen  in  Table  1 , the  ^ values  for  Ox  = a ■ 
and  the  values  for  <Jx  = -Oy  were  negligibly  small.  The  K-values  for  a 
45°  crack  and  the  K-values  for  a 0°  crack  are  plotted  in  Figure  12  . It  is 
seen  that  under  a one-to-one  biaxial  ratio, the  values  for  the  45°  crack 
agree  with  those  for  the  0°  crack.  On  the  other  hand,  the  absolute  ^ values 
for  the  45°  crack  agree  with  the  values  for  the  0°  crack  when  the  specimen 
is  subjected  to  a one-to-minus  one  biaxial  loading  condition.  As  for  the 
uniaxial  condition  (C^  = 0),  and  values  for  the  45°  crack  were 
approximately  equal  to  one-half  of  the  correspondingly  pure  Mode  1 and  pure 
Mode  2 stress  intensity  values.  The  expected  values,  according  to  the  angle 
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TABLE  1.  STRESS  INTENSITY  FACTORS  FOR  A 45°  CRACK  IN  A CRUCIFORM  SPECIMEN 
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Figure  11.  Finite  Element  Model  for  Angle  Crack 


crack  equations  given  in  [1],  are  indicated  by  the  dotted  lines  in  Figure  12. 
It  should  be  noted  that  the  minus  sign  for  K2  in  Table  1 was  only  a sign 
convention  used  to  indicate  the  slip  direction  of  a Mode  2 crack  (the  rela- 
tive motion  between  the  two  surfaces  of  a crack);  therefore,  only  the  abso- 
lute values  of  K2  are  plotted  in  Figure  12.  On  the  other  hand,  the  minus 
sign  for  in  Table  1 indicates  that  the  crack  was  not  open,  i.e.,  was 
actually  equal  to  zero. 


6.3  Crack  At  A Hole 


In  Section  6.1,  it  was  discussed  that  lateral  tension  or  compression 
stresses  would  not  affect  the  elastic  stress  intensity  factors  for  a central 
crack  in  a plate.  However,  if  a crack  (or  cracks)  is  coming  out  of  a circular 
hole  in  a plate,  compressive  loading  parallel  to  the  crack  can  cause  tensile 
Mode  1 stress  intensity  factors.  On  the  other  hand,  tensile  stresses 
parallel  to  the  crack  reduce  the  stress  concentrations  at  the  hole  and  thus 
reduce  the  crack  tip  stress  intensity  factor.  The  Bowie  solutions,  [22] 
have  provided  stress  intensity  factors  for  a single  crack  (or  cracks) 
emanating  frim  the  edge  of  a circular  hole,  in  an  infinitely  wide  plate, 
under  uniaxial  and  the  one-to-one  biaxial  loading  conditions.  For  any  other 
biaxial  load  ratios  (either  tension  combined  with  tension  or  tension  combined 
with  compression),  stress  intensity  factor  can  be  developed  by  using  super- 
positions of  the  uniaxial  and  biaxial  solutions  of  Bowie  as  illustrated  in 
Figure  13.  Following  this  superposition  logic,  stress  intensity  for  any 
biaxial  loading  combinations  can  be  expressed  as 


K = 


• CTx} 


B + 


(5) 


where  Bq  is  the  Bowie  factor  for  uniaxial  loading,  and  B^  is  the  Bowie 
factor  for  one-to-one  biaxial  loading,  and  a is  the  crack  length  on  either 
side  of  the  hole.  Letting  B =c^/°’y>  Equation  5 can  be  written  as 

K = OyjTfT  |(1  - B)  Bq  + B • B1  J (6) 

and  the  sign  for  B may  be  either  (+)  or  (-). 
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If  the  hole  is  inside  a finite  width  plate,  the  plate  dimensions 
influence  the  stress  concentration  at  the  edge  of  the  hole  and  thereby  varies 
the  stress  intensity  factor.  Referring  back  to  Section  6.1,  finite  width 
dimension  also  increases  the  crack  tip  stress  intensity  as  if  the  hole  was 
never  there.  Modifications  to  Equation  6 have  been  reported  in  the 
literature  (see  [20]  and  [23]  for  B=  0).  For  engineering  purposes, 

Equation  6 can  be  written  as 

K = | (1  - B)  Bq  + B • B1  j • F (7) 

For  a circular  hole  inside  a rectangular  sheet,  F can  be  a compounded 
factor  which  is  a product  of  the  Bowie  solution,  the  Howland  solution 
[24],  and  the  stress  intensity  factor  of  Equation  4 . The  validity  of 
these  assumptions  has  been  checked  out  by  experiment  and  by  cracked  finite 
element  analysis.  An  example,  taken  from  [25],  for  B ~ 0,  is  presented  in 
Figure  14.  There,  the  solid  line  is  the  computed  compound  factor.  The  test 
data  points  were  developed  by  using  the  backtracking  of  da/dN  data  technique 
as  described  in  [26].  Notice  that  the  experimental  data  points  at  a/r  >4.5 
(at  the  end  of  that  test)  were  significantly  higher  than  the  compounded 
factor  as  well  as  the  cracked  finite  element  analysis  results.  Realizing 
that  the  total  crack  length,  at  that  time,  was  larger  than  one  half  of  the 
specimen  width,  anti-plane  crack  tip  buckling  and/or  excessive  shear  lip 
might  have  developed  causing  the  alternation  of  the  crack  growth  behavior. 
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Figure  14.  Stress  Intensity  Factors  For  Two  Cracks  Coming 

Out  From  A Circular  Hole  In  A Finite  Width  Sheet 
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6.4  Curved  Crack 


It  has  been  demonstrated  by  elastic  analysis  [27]  [28]  that  a crack 

will  grow  straight  (stay  on  its  initial  path)  under  tension- compression 

biaxial  stresses,  but  the  crack  will  turn  away  from  its  initial  path  if  the 

biaxial  stress  ratio  is  larger  than  unity,  i.e.,  if  cr^  > CT  in  tension.  For 

a curved  crack  in  a biaxial  stress  field,  an  approximate  method  used  by 

Leevers,  et  al  [12]  can  be  used  to  compute  K.  Since  their  method  of 

analysis  primarily  deals  with  an  inclined  crack  (with  respect  to  either 

(7  or  a ),  it  would  be  necessary  to  compute  both  the  opening  mode  stress 
y x 

intensity  , and  the  sliding  mode  crack  tip  stress  intensity,  K^.  Therefore, 
their  equation  has  been  written  as 

K = F.  cr  JWa~+F1  <J  JWa~  (8) 

1 ly  y v y lx  x^  x 


K2  " F2y  <,y'S57+  F2x 


(9) 


where 

Flx 


0 <y\  . *2 n 

o-cos-j  Jsin 


a 7 a 5a  . 4a  , „ . 5a  . a . 2a ) 

^ ^cos  7 - cos  — - cos— sm  7 + 2sin7“-sirtjsin  7 J.  (10) 


i\  = 7 — — — icos  7 + cos  7^-  + cos-r~sin  — - 2sin7^sinysin^7! 

ly  4-cos?  cos^?  4 4 9 4 4 2 4 


(11) 


2x 


■3  a\  . \a 
J-cos-r  sm  , 


I . a .7a  . 5a  . 4a  „ . 5o  a . 2 a) 

<sm  7 - sm  7—  - sxn7-sm  7 - 2S1117— cos-jsm  7S 
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2y 


in  a 

J-cosy 

L 1 


|cos  J 


I . a . 7a  . 5ff  . 4o  5tt  a . 2a  I 

7 sin  7-  + sm  1-  sm-; — sm  y + isirrr — cos-rsm  — , 
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(13) 


and  the  definitions  for  a , a , and  01  are  given  in  Figure  15  where  a is  an 

x y 

angle  tangent  to  the  crack  at  an  instantaneous  point  on  the  crack  path;  i.e., 
the  angle  between  X-axis  and  the  line  connecting  two  consecutive  crack 
length  data  points. 
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When  Equations  8 and  9 are  used  to  correlate  crack  growth  rate  and/or 
residual  strength  test  data,  it  is  necessary  to  adopt  a failure  criterion 
(to  an  equivalent  K value)  accounting  for  the  combined  effects  of  and 
at  the  crack  tip.  There  are  numerous  failure  criteria  available  in  the 
literature,  e.g.,  [29]  to  [31].  In  this  study,  the  following  possibilities 
were  considered: 

K = (K^  4-  (14) 

based  on  Irwin's  theory  of  fracture  [29],  and 


K * K2  + K- 


(15) 


derived  from  experimental  data  on  fracture  testing  of  aluminum  alloys  and 
4340  steel  [16  and  32]. 


A simplified  method  for  calculating  K was  also  considered.  This  method 

simply  calculates  K values  in  the  following  manner.  It  has  been  determined 

by  cracked  finite  element  analysis  that  the  K values  at  all  crack  lengths 

in  the  flat  circle  of  the  cruciform  specimen  (for  o <a  ) are  essentially 

x y 

the  same  as  for  an  infinite  sheet  under  uniaxial  tension.  Therefore, 

K = a y va.  Here,  o is  taken  to  be  the  maximum  stress  perpendicular  to  the 
crack  at  each  small  segment  on  the  crack  growth  path, as  schematically 
illustrated  in  Figure  16.  It  is  seen  that  o changes  as  the  crack  propagates. 
It  equals  a at  the  beginning  of  the  test  and  gradually  changes  to  when 
the  crack  has  completely  turned  90°.  The  crack  length  "a"  used  in  this 
method  of  calculation  is  rather  arbitrary.  As  indicated  in  Figure  16,  the 
distance  connecting  the  center  of  the  specimen  and  middle  point  of  a da 
segment  is  used  as  an  equivalent  "a".  In  this  approach,  since  an  effective 
stress  which  opens  the  crack  is  used  for  computing  an  effective  K,  it  may  be 
considered  that  the  term  has  been  lumped  into  the  calculation. 
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SECTION  VII 


CRACK  TIP  PLASTIC  ZONE  SIZES 


A crack  tip  plastic  zone  is  formed  during  the  upward  excursion  of  each 
loading  cycle.  After  removal  of  the  applied  load,  or  in  the  valley  of  a 
loading  cycle,  there  is  a residual  plastic  zone  remaining  at  the  crack  tip. 

The  sizes  of  these  crack  tip  plastic  zones  are  functions  of  crack  length, 
stress  level,  biaxial  stress  ratio  and  the  state  of  stress  triaxiality  at  the 
crack  tip  (e.g.,  plane  stress  or  plane  strain).  In  a constant  amplitude 
loading  case,  the  crack  growth  rate  behavior  may  be  related  to  the  crack  tip 
plastic  zone  size  at  the  peak  of  a load  cycle.  However,  in  a variable 
amplitude  loading  situation,  taking  a high-low  block  loading  sequence,  for 
example,  it  is  commonly  believed  that  the  residual  plastic  zone  produced  by 
the  high  load  is  the  source  affecting  the  crack  growth  retardation  behavior. 

In  this  study,  the  "loaded"  crack  tip  plastic  zones  were  determined  by 
elastic-plastic  finite  element  (NASTRAN)  analyses,  and  also  by  taking 
interferometry  photographs  at  the  crack  tip  during  a cyclic  crack  growth  test. 
The  "residual"  crack  tip  plastic  zones  were  determined  by  interferometry 
measurements  and  also  by  backtracking  of  crack  growth  data  obtained  from 
periodically  overloaded  specimens.  The  details  about  these  techniques  and 
the  results  obtained  from  each  of  these  techniques  are  discussed  below. 

7 . 1 Finite  Element  Analysis 

Elastic-plastic  NASTRAN  analyses  were  conducted  to  determine  crack  tip 
plastic  zone  sizes  in  a biaxially  loaded  cruciform  specimen.  The  finite 
element  model  shown  in  Figure  17  represents  a quarter  of  a circular  section 
of  3.5  inches  in  radius.  The  location  of  the  nodal  points  on  the  outer 
circumference  of  this  model  matches  those  shown  in  Figure  6.  The  arrangement 
and  size  of  the  elements  (not  cracked)  along  the  crack  line,  especially  in 
the  vicinity  of  the  crack  tip,  were  modified  for  each  crack  size. 
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For  example,  fewer  plastic  elements  and  finer  elements  were  used  around  the 
crack  tip  for  shorter  crack  lengths  because  smaller  crack  tip  plastic  zones 
had  been  expected,  and  vice  versa.  The  elastic-plastic  boundary  arrangement 
such  as  those  shown  in  Figure  17  was  used  for  the  cases  of  a 1.5  inch  crack 
in  a 7075-T7351  specimen  and  a 1.0  inch  crack  in  a 2024-T351  specimen. 

The  general  procedure  for  these  analyses  was  to  apply  the  required  loads 
(for  a desired  stress  level  and  biaxial  ratio)  onto  the  old  model  (coarse 
grid,  all  elastic,  see  Figures  6A  and  6B ) without  a crack.  The  loads  at  the 
nodal  points  around  the  3.5  inches  radius  were  determined  and  these  loads 
were  then  applied  onto  the  new  model  (Figure  17)  with  cracks.  The  same 
loads  could  be  used  for  different  crack  lengths  as  long  as  their  biaxial 
stress  ratio  and  magnitude  were  the  same  thereby  saving  considerable  computer 
time . 


Taking  the  stresses  from  the  NASTRAN  print  out,  the  Von  Mises  stress  in 
each  plastic  element  was  computed.  An  element  was  considered  yielded  when 
its  Von  Mises  stress  exceeded  the  elastic  proportional  limit  of  the  material 
(58  ksi  for  7075-T7351  and  52  ksi  for  2024-T351).  A plastic  zone  Was  defined 
by  gathering  a group  of  the  yielded  elements.  Plastic  zone  contour  maps  for 
five  loading  cases  (O^/cr^  = 0,  441.5,  +1.0)  at  30  ksi  for  both  7075-T7351  and 
2024-T351  materials,  at  seven  crack  lengths  (a  = 0.25,  0.5,  1.0  and  1.5  for 
7075-T7351  and  a = 0.25,  0.5  and  1.0  for  2024-T351  specimens)  are  presented 
in  Figures  18  through  24.  Note  the  dotted  circle  in  Figure  18  which  maps  a 
very  small  plastic  zone  is  only  an  estimated  value.  In  these  cases  the 
elements  directly  surrounding  the  crack  tip  were  not  yielded  (based  on  the 
calculated  Von  Mises  stress),  and  the  plastic  zone  was  assumed  to  be  one-half 
of  the  element  sizes  since  the  Von  Mises  stresses  computed  were  based  on  the 
stresses  at  the  centroid  of  each  element.  As  for  the  huge  plastic  zones  in 


Figures  21  and  24,  the  incompleted  zones  for  cr^ 


CT  were  due  to  the  fact 

y 


that  these  plastic  zones  were  so  big  and  their  actual  sizes  were  beyond  the 


boundary  of  the  elastic  and  plastic  elements  (see  Figure  17). 
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Crack  Tip  Plastic  Zone  for  2024-T351  at  30  ksi  (NASTRAN  analysis) 


Crack  tip  plastic  zone  sizes  at  lower  stress  levels  were  also  determined 
by  elastic-plastic  NASTRAN  analyses.  The  plastic  zone  dimensions  for  all  the 
cases  were  measured  from  each  contour  map  (reduced  from  the  NASTRAN  data)  and 
are  listed  in  Tables  2 and  3.  The  methods  used  for  the  measurements  are 
described  in  the  following. 


The  shape  and  the  dimension  of  a crack  tip  plastic  zone 
crack  length,  stress  level,  and  the  mechanical  properties  of 
For  a crack  under  uniaxial  tension,  in  a plane  stress  state, 
equations  of  Irwin,  and  Rice  [33],  are  respectively 


1 2 
r = - (K  /F  ) 
p 7 t max  ty 


are  functions  of 
a given  material, 
the  classical 


(16) 


and 


r 

P 

However,  the  characteristic  d 
seen  in  Figures  18  to  24 , the 


7T 


= o (K  /F 
8 max 


(17) 


imension  r has  never  been  clearly  defined.  As 
P 

plastic  zone  shapes  varied  among  all  the  cases. 


The  commonly  familiar  plastic  zone  shape  (the  butterfly  shape)  is 
schematically  shown  in  Figure  25.  Here  we  define  the  dimension  r to  be  the 
largest  distance  between  the  crack  tip  and  the  border  of  the  plastic  zone. 
Then  we  compute  an  equivalent  dimension  r representing  the  diameter  of  a 
circle.  The  area  of  such  an  imaginary  circle  is  set  equal  to  the  total  area 
of  the  butterfly. 
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TABLE  2.  NASTRAN  CRACK  TIP  PLASTIC  ZONE  FOR  7075-T7351  ALLOY 
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TABLE  3.  NASTRAN  CRACK  TIP  PLASTIC  ZONE  FOR  2024-T351  ALLOY 


a 

(Inch) 

0 

y 

(Ksi) 

r 

P 

(Inch) 

r 

P 

(Inch) 

/ K \2 

1 max  1 

\Fty  / 

(Inch) 

ax 

la 

y 

.25 

30.0 

.11 

.172 

.283 

0 

.5 

30.0 

.19 

.244 

.565 

- .5 

24.0 

.1625 

.19 

.362 

1.0 

30.0 

.413 

.562 

1.131 

.25 

30.0 

.0742 

.108 

.283 

.5 

.5 

20.0 

.1 

.108 

.251 

.5 

30.0 

.141 

.194 

.565 

1.0 

30.0 

.306 

.45 

1.131 

1.0 

20.0 

.169 

.202 

.503 

.25 

30.0 

.109 

.144 

.283 

1.0 

.5 

30.0 

.13 

.216 

.565 

1-0 

24.0 

.269 

.316 

.724 

1.0 

30.0 

.36 

.45 

1.131 

.25 

30.0 

.141 

.24 

.283 

— 

.5 

.5 

30.0 

.277 

.428 

.565 

.5 

20.0 

.171 

.18 

.251 

1.0 

20.0 

.313 

.36 

.503 

1.0 

30.0 

.918 

1.012 

1.131 

.25 

30.0 

.453 

.594 

.283 

-1 

.0 

.5 

30.0 

1.6 

1.756 

.565 

.5 

20.0 

.33 

.561 

.251 
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Both  dimensions,  and  r , are  listed  in  Tables  2 and  3 , and  are 

plotted  in  Figures  26  to  29  as  functions  of  a normalized  plasticity  density 

2 

parameter,  ^max/Fty)  • Also  in  these  figures  are  the  plastic  zone  sizes 
computed  by  using  Equations  16  and  17  . These  theoretical  values  are 
included  here  to  compare  with  the  NASTRAN  results  for  the  uniaxial  condition 
(cr  =0).  It  also  provides  some  indication  of  the  relative  sizes  of  the 

X 

crack  tip  plastic  zones  for  cracks  subjected  to  other  biaxial  loading 
conditions.  Facts  that  can  be  observed  from  Figures  26  to  29  are  listed  in 
the  following; 


1.  For  the  same  data  point  the  r dimension  is  usually,  although  not 
always,  slightly  larger  than  the  r^  dimension. 

2.  The  plastic  zone  sizes  for  biaxial  ratios  of  0.5  and  1.0  are 
approximately  the  same  and  are  insignificantly  smaller  than  those 
for  the  uniaxial  tension. 

3.  The  plastic  zone  sizes  for  tension-compression  biaxial  ratios  are 
significantly  larger  than  those  for  the  uniaxial  tension  and  tension- 
tension  biaxial  conditions;  the  higher  the  tension-compression  ratio, 
the  larger  plastic  zone  size. 

4.  The  r dimensions  for  the  uniaxial  tension  cases  fit  in  between  the 

P 

theoretical  values  of  Rice  and  Irwin  whereas  the  r dimensions  are 

P 

slightly  larger  than  the  Rice  plastic  zones. 

5.  For  the  2024-T351  specimen,  the  relationship  between  r^  (or  r^)  and 
the  plasticity  density  parameter,  (Kmax/F  )^,  is  linear  for  at 
least  four  out  of  the  five  biaxial  loading  conditions  (excepting 

B = -1.0  which  is  questionable).  However,  for  the  7075-T7351 

specimen  at  a given  biaxial  ratio,  this  relationship  is  linear  up  to 

2 

a certain  (K  /F.  ) value.  For  the  -1.0  biaxial  ratio  cases,  the 
max  ty  ’ 

crack  tip  plastic  zone  for  the  2024-T351  specimen  is  always  larger 
than  those  for  the  7075-T7351  specimen.  As  for  the  other  four 

biaxial  ratios,  the  plastic  zone  sizes  for  both  2024-T351  and 

_ 2 

7075-T7351  are  the  same  as  long  as  the  r (or  r ) and  (K  /F.  ) 

6 p p max  ty 

relationships  for  both  specimens  are  linear. 
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Figure  26.  Crack  Tip  Plastic  Zone  for  7075-T7351  (NASTRAN  analysis) 
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Figure  27.  Crack  Tio  Plastic  Zone  for  7075-T7351  (NASTRAN  analysis) 
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Figure  28.  Crack  Tip  Plastic  Zone  for  2024-T351  (NASTRAN  analysis) 
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Figure  29.  Crack  Tip  Plastic  Zone  for  2024-T351  (NASTRAN  analysis) 


Actually,  plastic  zone  size  nonlinearity  exists  in  both  7075-T7351 

and  2024-T351  specimens.  The  reason  for  its  not  being  observed  in  the 

2024-T351  specimen  was  due  to  the  fact  that  the  largest  crack  size  in  the 

2024-T351  specimen  was  one-inch  (2a  = 2 inches),  Whereas  the  largest  crack 

size  in  the  7075-T7351  specimen  was  1.5  inches  (2a  = 3 inches).  For  a 

given  stress  level,  the  K /F,_  ratios  for  these  two  cases  (a  2-inch  crack  in 
b max  ty 

a 2024-T351  specimen  and  a 3-inch  crack  in  a 7075-T7351  specimen)  are 

approximately  the  same.  Taking  a 30  ksi  applied  stress  level  for  example, 

2 

the  (Kmax/F ty ) value  is  approximately  equal  to  1.15  and  the  corresponding 
plastic  zone  size  for  the  7075-T7351  specimen  is  significantly  larger  than 
those  for  the  2024-T351  specimen.  The  significantly  larger  crack  dimension 
in  the  7075-T7351  specimen  (3  inches  long  as  opposed  to  a 2 inch  crack  in 
the  2024-T351  specimen)  might  have  induced  some  geometric  effect  on  the 
development  of  a crack  tip  plastic  zone. 


7.2  Interferometry  Measurement 

For  the  interferometry  tests,  photoimprinted  grids  were  not  used  as  in 
other  tests.  Instead  a transparent  grid  (of  identical  design  to  that  used 
for  photoimprinting)  was  mounted  in  a fixture  attached  to  the  specimen. 

This  fixture  was  so  designed  to  ensure  repeatable  positioning  (for  each 
interferometer  scan,  this  was  necessarily  removed). 

A Wild  M20  interferometer  was  adapted  for  use  in  the  biaxial  frame. 

The  interferometer  base  was  positioned  on  the  top  cover  plate  of  the  test 
frame  with  the  objective  lens  protruding  down  through  the  access  hole  of  the 
cover  plate.  For  each  plastic  zone  mapped,  the  interferometer  was  positioned 
parallel  to  the  crack  direction  and  subsequent  interference  photographs  taken. 
The  Wild  M20  was  limited  to  a minimum  magnification  of  45x.  This  required 
that  the  plastic  zone  be  mapped  by  taking  several  photographs  which  were 
later  assembled  into  a composite. 
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Some  clear  limitations  on  the  use  of  this  interference  microscope 
become  apparent. 

1.  A high  degree  of  subjectivity  was  involved  in  setting  up  the 
equipment,  i.e.,  there  was  no  absolute  criterion  for  the  adjust- 
ment required  to  generate  the  interference  patterns.  This  was 
particularly  true  for  large  plastic  zone  where  no  "flat"  surface 
can  be  discerned  at  45 x.  The  actual  adjustments  made  were  often 
done  on  a trail/error  basis,  i.e.,  the  entire  zone  was  mapped  and 
remapped  until  a coherent  composite  was  attained. 

2.  For  the  plastic  zone  measured  at  load,  the  vibrations  from  the  test 
frame  made  it  extremely  tedious  to  produce  decent  consistent 
photographs . 

Interferometry  photographs  were  taken  from  30  specimens.  In  each 
specimen,  photographs  were  taken  at  the  crack  tip  at  two  to  four  crack 
dimensions  ranging  from  0.2  to  2.0  inches  of  total  crack  length.  Whenever 
possible,  photographs  were  taken  at  zero  load  (after  unloaded  from  maximum 
load)  and  then  at  peak  load  (after  reloaded  to  maximum  load  from  zero 
load).  The  first  picture  taken  represented  the  residual  plastic  zone.  For 
those  periodically  overloaded  specimens,  the  interferometry  photographs 
were  actually  taken  at  zero  load  prior  to,  and  after,  the  overload.  The 
first  picture  taken  represented  the  residual  plastic  zone  produced  by  the 
constant  amplitude  stress  and  the  second  picture  represented  the  residual 
plastic  zone  of  the  overload. 

Three  patterns  of  interferometry  fringes  are  shown  in  Figures  30  to 
32.  These  fringe  patterns  are  typical  of  all  the  photographs  taken  from 
all  the  specimens.  Figure  33  is  a schematic  representation  of  these  fringe 
patterns.  The  pattern  of  Figure  30,  or  Figure  33(a),  is  usually  associated 
with  constant  amplitude  at  a lower  stress  intensity  level  (either  "loaded" 
or  "unloaded")  whereas  the  other  two  patterns  usually  appeared  after  an 
overload.  For  the  constant  amplitude  at  higher  K level  cases,  the 
interferometry  fringe  pattern  was  somewhat  in  between  Figures  30  and  31. 
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Each  of  the  interferometry  contours  actually  measures  the  depth  of  the 
surface  depression  of  the  metal  sheet  around  the  crack,  i.e.,  each  fringe 
represents  a certain  amount  of  vertical  displacement.  When  the  contour 
line(s)  becomes  out  of  focus,  or  diverges,  the  material  outside  of  that 
contour  is  probably  not  deformed.  Therefore,  it  can  logically  be  said  that 
the  last  visible  contour  (e.g.,  line  A in  Figures  30  and  32)  divides  the 
plastically  deformed  and  undeformed  materials  and  thereby  maps  out  the 
crack  tip  plastic  zone.  As  for  those  fringe  patterns  shown  in  Figure  31 
or  Figure  33(b),  the  division  for  a plastic  zone  is  not  clear.  Comparing 
Figures  33(b)  and  33(c)  the  main  difference  between  these  two  patterns  is 
that  all  the  contour  lines  in  Figure  32  or  Figure  33(c),  immerse  into  the 
crack  tip  whereas  many  of  the  outer  contour  lines  in  Figure  31  or 
Figure  33(b),  do  not  close  all  the  way  to  the  crack  tip.  Therefore,  some 
contour  lines  which  start  to  turn  around  and  become  parallel  to  the  crack 
line  were  selected  to  represent  the  crack  tip  plastic  zone.  When  uncertainty 
arose,  two  contour  lines  were  selected,  one  representing  a minimum  r^  value 
and  another  representing  a possible  r^  value  (see  Figure  31). 

The  interferometry  photographs  taken  from  all  30  specimens  are  compiled 
in  Volume  III  of  this  report.  The  crack  tip  plastic  zone  sizes  (the  r^ 
values)  reduced  from  each  of  these  photographs  are  presented  in  Figures  34 
to  48  comparing  the  "loaded"  and  the  "unloaded"  plastic  zone  sizes  at  each 
biaxial  stress  ratio.  Comparison  of  the  analytically  (elastic-plastic 
NASTRAN)  and  the  experimentally  (interferometry)  determined  plastic  zone 
sizes  (at  load)  are  presented  in  Figures  49  to  55.  It  can  be  observed 
from  these  figures  that: 

1.  Excepting  the  one-to-minus  one  biaxial  ratio  case,  the 

experimentally  determined  crack  tip  plastic  zone  (at  peak  load) 
correlates  with  those  determined  by  using  elastic-plastic  finite 
element  analysis  technique. 
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7075-T7351,  a = 0.5  inch 

cr  = 0,  a =12  ksi,  unloaded 
x y 

Magnification:  73.0  X 


Figure  31.  Interferometry  Photograph  of  Cruciform  Specimen  No.  7-114,  TC  No.  99 
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2 0 2 <4 - T 7 3 5 1 , a = 0.26  inch 
c = 0,  o =20  ksi  over- 
x y 

loaded.  Magnification:  70.3  X 


Figure  32.  Interferometry  Photograph  of  Cruciform  Specimen  No.  2-24,  TC  No.  101 


58 


2. 


Excepting  the  one-to-one  biaxial  case,  the  unloaded  plastic  zone 

sizes  for  all  other  biaxial  ratios  are  generally  scattered  around 

the  Irvin/Rice  plastic  zone  equations.  However,  when 
2 

(Kmax^ty)  ^ 0.4  the  interferometry  plastic  zones  are  significantly 
and  consistently,  in  all  the  cases,  smaller  than  the  Irwin/Rice 
plastic  zones.  This  might  be  an  indication  of  the  limitation  of 
the  interferometry  technique. 


3.  For  the  one-to-one  biaxial  ratio,  both  the  loaded  and  unloaded 

2 

plastic  zones,  at  any  (K  / F.  ) value,  are  smaller  than  the 

K J max  ty 

Irwin/Rice  plastic  zones. 


4.  The  plastic  zones  at  load  are  generally  bigger  than  the  unloaded 
plastic  zones.  However,  the  exact  ratios  between  these  two  types 
of  plastic  zones  cannot  be  determined  at  this  time. 


7 . 3 Backtracked  Residual  Plastic  Zones 

The  existing  crack  growth  retardation  models  generally  assure  that  all 
the  crack  growth  interaction  activities  will  take  place  inside  a load 
interaction  zone  (equivalent  to  a crack  tip  plastic  zone)  created  by  the 
high  load.  Once  the  crack  has  propagated  through  this  interaction  zone,  a 
normal  crack  growth  rate  will  be  resumed.  Therefore,  it  will  be  possible 
to  trace  the  extent  (the  size)  of  such  an  interaction  zone  from  actual 
crack  growth  rate  data.  Taking  the  crack  growth  rate  data  for  the  period- 
ically loaded  specimens  (see  Figures  110  to  118  in  Section  9.5.1)  for 
example,  in  part  (a)  of  each  figure,  all  da/dN  data  points  including  those 
affected  by  load  interactions  are  plotted  as  a function  of  an  apparent 
(uncorrected)  K level.  Since  each  retarded  da/dN  data  point  associates 
with  a crack  length,  according  to  the  hypothesis  stated  earlier,  it  would 
be  logical  to  assume  that  the  size  of  the  load  interaction  zone  is  equal  to 
the  difference  between  the  crack  length  at  overload  and  the  effective  crack 
length  associated  with  the  last  traceable  retarded  da/dN  data  point. 


60 


2024-T351 
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Figure  34.  Interferometry  Measurement  of  Crack  Tip  Plastic  Zone  Sizes 
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Figure  35.  Interferometry  Measurement  of  Crack  Tip  Plastic  Zone  Sizes 
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Figure  36.  Interferometry  Measurement  of  Crack  Tip  Plastic  Zone  Sizes 


2024-T351 
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Figure  37.  Interferometry  Measurement  of  Crack  Tip  Plastic  Zone  Si 
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Figure  38.  Interferometry  Measurement  of  Crack  Tip  Plastic  Zone  Sizes 
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Figure  39.  Interferometry  Measurement  of  Crack  Tip  Plastic  Zone  Si 


2024-T351,  o = -0.35  a 
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Figure  40.  Interferometry  Measurement  of  Crack  Tip  Plastic  Zone  Sizes 
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Figure  41.  Interferometry  Measurement  of  Crack  Tip  Plastic  Zone  Sizes 
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Figure  42.  Interferometry  Measurement  of  Crack  Tip  Plastic  Zone  Sizes 
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Figure  43.  Interferometry  Measurement  of  Crack  Tip  Plastic  Zone  Sizes 
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Figure  44.  Interferometry  Measurement  of  Crack  Tip  Plastic  Zone  Sizes 
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Figure  46.  Interferometry  Measurement  of  Crack  Tip  Plastic  Zone  Sizes 
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Figure  47.  Interferometry  Measurement  of  Crack  Tip  Plastic  Zone  Sizes 
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Figure  48.  Interferometry  Measurement' of  Crack  Tip  Plastic  Zone  Sizes 
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Figure  49.  Comparison  of  Interferometry  and  NASTRAN  Crack  Tip  Plastic 
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Figure  50.  Comparison  of  Interferometry  and  NASTRAN  Crack  Tip  Plastic 
Zone  Sizes 


T351 


Figure  52.  Comparison  of  Interferometry  and  NASTRAN  Crack  Tip  Plastic 
Zone  Sizes 
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Figure  54.  Comparison  of  Interferometry  and  NASTRAN  Crack  Tip  Plastic 
Zone  Sizes 


82 


Figure  55.  Comparison  of  Interferometry  and  NASTRAN  Crack  Tip  Plastic 
Zone  Sizes 


The  backtracked  plastic  zone  sizes  and  the  residual  plastic  zones 
(determined  by  interferometry  measurement  after  each  overload)  are  plotted 
in  Figures  56  to  61  for  comparison.  Note  in  these  figures,  the  Kmax  is 
computed  using  the  overload  stress  and  the  crack  length  at  overload.  Since 
these  Tp  values  were  obtained  from  the  interpretation  of  interferometry 
fringe  patterns  or  backtracking  of  da/dN  values,  when  uncertainity  arose 
(note  the  uncertainties  came  from  the  appearance  of  the  fringe  pattern  or 
the  scatter  of  the  material  da/dN  data),  a minimum  Tp  value  and  a possible 
Tp  value  were  determined  and  both  values  are  plotted  in  the  figure.  Also 
due  to  the  fact  that  in  some  cases  the  crack  had  grown  longer  on  one  side 
and  shorter  on  the  other  side,  the  plastic  zone  sizes  for  the  left  crack  tip 
and  the  right  crack  tip  were  treated  as  two  independent  data  points.  The 
interferometry  photographs  were  always  taken  from  the  right  side  of  the 
crack. 


Examination  of  Figures  56  to  61  reveals  the  following: 

1.  At  lower  K levels,  both  the  interferometry  residual  plastic  zones 
and  the  backtracked  residual  plastic  zones  are  close  to  the  sizes 
computed  by  Equations  16  and  17,  for  monotonic  loading  condition. 

2.  At  higher  K levels,  the  backtracked  residual  plastic  zones  are 
significantly  larger  than  those  calculated  using  the  equations  of 
Irwin  and  Rice. 

3.  At  higher  K levels  the  backtracked  residual  plastic  zones  for  the 
2024-T351  specimens  seem  to  be  larger  than  those  for  the 
7075-T7351  specimens. 

4.  The  effect  of  biaxial  loading  ratios  on  residual  plastic  zone 
sizes  cannot  be  identified  (note  the  biaxial  ratios  for  these 
specimens  were  0,  +0.5).  However,  it  seems  that  the  Tp  values 
for  the  -0.5  biaxial  ratio  cases  (in  either  one  of  the  two 
aluminum  alloys)  could  possibly  be  larger  than  those  in  the  other 
two  biaxial  conditions;  i.e.,  it  might  be  an  indication  of  a 
biaxial  loading  effect. 
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5.  As  indicated  by  the  backtracked  plastic  zones,  the  residuaal 

plastic  zones  increased  rapidly  when  the  plasticity  density 

2 

parameters,  IKmax/F^  ) , exceed  0.4;  this  is  a phenomenon  much 
similar  to  those  previously  observed  in  the  "loaded"  plastic 
zones  (determined  by  NASTRAN) . At  this  moment,  it  is  not  certain 
whether  these  nonlinearities  were  indications  of  a real  physical 
phenomenon  or  just  special  cases  that  had  been  associated  with 
the  cruciform  specimen  configuration. 
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S-T7351,  TEST  CASE  99, 
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Figure  56.  Experimentally  Determined  Residual  Crack  Tip  Plastic  Zone 


7075-T7351,  TEST  CASES  103  AND  104, 
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Figure  57.  Experimentally  Determined  Residual  Crack  Tip  Plastic  Zone 
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Figure  58.  Experimentally  Determined  Residual  Crack  Tip  Plastic  Zone 


2024-T351,  TEST  CASES  101  AND  102, 


o 


88 


Figure  59.  Experimentally  Determined  Residual  Crack  Tip  Plastic  Zone 


2024-T351  ,TEST  CASES  108  AND  1 
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Figure  60.  Experimentally  Determined  Residual  Crack  Tip  Plastic  Zone 


2024-T351 , TEST  CASE  110, 
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Figure  61.  Experimentally  Determined  Residual  Crack  Tip  Plastic  Zone 


SECTION  VIII 


CRACK  GROWTH  ANALYSIS  METHODOLOGY 


In  this  section,  the  mechanics  of  crack  growth  under  constant  amplitude 
and  variable  amplitude  loading  conditions  are  discussed.  Analytic  procedures 
for  calculating  the  predicted  crack  growth  histories  are  outlined.  Three 
computer  programs  were  constructed.  The  first  program  dealt  with  the  crack 
growth  rate  and  the  crack  growth  direction  for  a crack  propagating  in  a 
biaxial  stress  field.  The  second  program  adopted  the  Modified  Willenborg 
model  to  compute  crack  growth  history  for  a crack  under  variable  amplitude 
loading  conditions. 

Taking  advantage  of  the  actual  crack  growth  rate  variations  (for  a 
crack  inside  an  overloaded  zone)  and  the  "crack  closure"  phenomenon  reported 
by  Von  Euw  [34]  and  Elber  [35],  proper  assumptions  were  made  to  tie  these 
two  crack  propagation  mechanisms  together  into  one  unified  mathematical 
procedure.  This  new  crack  growth  retardation  model  was  incorporated  into 
the  third  computer  program. 

All  the  methodologies  discussed  in  the  following  sections  are  correlated 
with  the  experimental  data  presented  in  Section  9. 


8.1  Fatigue  Crack  Growth  Under  Constant  Amplitude  Cyclic  Loading 

As  will  be  shown  in  the  later  sections,  for  negative  biaxial  ratios, 
and  for  positive  biaxial  ratios  up  to  1.0,  the  crack  will  grow  straight 
(stay  on  its  initial  path)  and  crack  growth  rates  are  not  affected  by  the 
lateral  stresses  (parallel  to  the  crack).  Therefore,  crack  growth  analysis 
can  be  done  by  using  conventional  methods  as  long  as  the  K factors  (in  terms 
of  a ) are  properly  determined. 
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For  biaxial  ratios  larger  than  1.0,  available  analytical  results  and 
the  present  experimental  results  show  that  the  crack  will  turn  away  from  its 
initial  plane.  There  are  two  problems  in  conducting  a computerized  crack 
growth  path  analysis.  The  first  problem  is  how  to  determine  the  K and  the 
da/dN  values,  and  the  second  problem  is  how  to  predict  the  subsequent  crack 
growth  direction  (from  a current  crack  tip  location).  Stress  intensity 
factors  for  a curved  crack  have  been  discussed  in  Section  6.4.  Based  on  the 
test  results  of  all  the  tests  conducted  on  negative  biaxial  ratios,  and 
positive  biaxial  ratios  up  to  1.0,  it  can  be  assumed  that  crack  growth  rate 
will  be  the  same  for  tension-high  tension  biaxial  ratios.  Therefore,  the 
criterion  for  determining  the  accuracy  of  a correlation  is  that  there  should 
be  no  effect  on  crack  growth  rate  if  K is  correctly  determined.  Consequently, 
whether  the  K is  determined  by  the  vector  sum  of  and  K„  (Equation  14), 
or  by  just  adding  and  K-  (Equation  15),  or  by  using  the  "equivalent  a" 
method  (Figure  16);  the  best  correlation  of  all  should  show  a da/dN  versus  K 
curve  identical  to  the  material  baseline  crack  growth  rate  curve  (or  any  crack 
growth  rate  curve  for  B -^>1.0).  As  will  be  discussed  in  Section  9.2.2, 
Equation  15  is  the  best  candidate  and  the  "equivalent  a"  method  is  the 
second  best  candidate  for  determining  K values  along  a curved  crack.  At  this 
point,  one  might  ask  if  should  vanish  after  the  crack  starts  to 
propagate.  The  following  discussions  may  clear  up  the  pros  and  the  cons  in 
the  mechanics  of  crack  growth  under  mixed-mode  conditions. 


For  a stationary  inclined  crack  in  uniaxial,  as  well  as  biaxial, 

stress  field,  both  and  exist  at  the  crack  tip.  The  classical  theory  of 

fracture  mechanics  [30]  has  predicted  that  the  crack  will  rotate  immediately 

on  propagation  and  propagates  in  the  direction  where  is  maximum.  The  angle 

of  initial  crack  extension  is  negative  (clockwise)  as  referenced  to  the 

original  crack  line.  Numerous  investigations,  e.g.,  [36]  on  uniaxial  tension 

and  [37]  on  biaxial  tension-compression,  have  shown  that  the  magnitude  of  ^ 

decreases  rapidly  after  a very  small  amount  of  crack  extension.  Crack  tip 

stress  intensity  at  various  points  in  the  remaining  portion  of  the  crack  path 

will  be  predominently  K^.  However,  the  classical  theory  is  only  applicable 

to  loading  conditions  such  that  the  magnitude  of  cr  is  always  lower  than  cr  , 

x y 

i.e.,  cr  is  in  compression,  cr  = o,  or  cr  ^ a • 
x 1 x x ^ y 
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When  a crack  is  subjected  to  a one-to-one  biaxial  ratio,  the  stresses  in  all 
directions  are  the  same;  therefore,  the  crack  will  grow  straight  regardless 
of  its  original  orientation.  If  ax  is  larger  than  cry  in  tension,  recent 
analytical  results  [28]  have  shown  that  the  crack  will  turn  toward  the 
Y-axis.  Regardless  of  whether  the  original  crack  plane  is  either  parallel 
or  inclined  to  CX^,  the  angles  of  initial  crack  extension  will  be  positive 
(counterclockwise)  and  the  crack  will  end  up  normal  to  0^,  the  most 
dominating  stresses  where  maximum  is  enhanced.  Analytic  solutions 
relating  the  angle  of  initial  crack  extension  and  the  original  orientation 
of  the  crack  is  available [28] . However,  there  has  been  no  analysis 
concerning  the  status  of  stress  intensity  variations  after  the  crack  has 
extended.  Referring  back  to  Equations  8 and  9,  since  Leevers  treated  the 
S-shaped  crack  geometry  as  being  a simulation  of  Sih's  arch  crack  problem 
[38] , the  general  implication  is  that  as  long  as  the  crack  maintains  an 
equivalently  arched  shape,  both  and  K2  exist  at  the  crack  tip.  Further- 
more, since  the  crack  propagation  angle  is  positive  in  all  the  cr  > or 

x y 

loading  cases,  the  curved  crack  will  remain  as  S-shaped  after  each  increment 
of  crack  growth.  Therefore,  it  is  anticipated  that  both  and  will 
appear  at  the  beginning  of  each  loading  cycle. 

A schematic  illustration  elucidating  the  major  steps  in  the  computerized 
routine  for  growing  a curved  crack  is  given  in  Figure  62.  Suppose  the  crack 
is  originally  set  at  »o  degrees  in  respect  to  the  X-axis.  Making  use  of  the 
crack  angle  and  biaxial  ratio  relationships  given  by  Eftis  and  Subramonian 

[28] , i ,e , , 

1.75  j D1*  (3Cos  6 -1)  - 3D2Sine|  + 4D3Sirt|cos  6 = 0 (18) 

With  Dx  = -3(l-B)Sina -cos(90°-o) 

D2  = B + (l-B)Sin2(90°-o) 

D3  = (1-B)Cos2(90°-o) 

one  may  obtain  a crack  propagating  angle,  6q,  which  defines  the  direction 
that  the  crack  will  takeoff.  At  the  same  time,  a parallel  calculation  has 
to  be  made  determining  the  K value  at  the  crack  tip  (Equation  15  was  used). 
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Figure  62.  Definitions  for  Effective  Crack  Angle  (a) 
and  Crack  Takeoff  Angle  (8)  in  Propagation 
of  a Curved  Crack. 
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Subsequently,  the  "da"  value  for  that  load  cycle  can  be  determined  through 

a crack  growth  rate  equation  which  presents  the  baseline  crack  growth  rate 

behavior.  Since  0Q  is  known,  the  dax  and  da^.  components  can  be  computed 

and  added  to  the  previous  ax  and  a^.  values,  thereby  determining  the  new 

crack  tip  location.  The  next  step  of  crack  extension  can  be  predicted  by 

repeating  the  whole  routine  of  calculations  (i.e.,  determining  the  9^  and 

the  new  da  and  da  values).  It  should  be  noted  that  the  Eftis  criterion, 
x y 

Equation  18,  was  originally  solved  for  a straight  crack,  not  a kinked  (or 
curved)  crack.  The  effectiveness  of  Equation  18,  for  predicting  the  new 
9' s,  depends  on  how  to  choose  a new  a after  each  step  of  crack  extension. 

In  the  present  computer  program,  all  the  new  a's  are  determined  by  using  a 
fixed  origin  system;  i.e.,  an  effective  new  a is  defined  to  be  the  angle 
between  the  X-axis  and  the  line  connecting  the  origin  and  the  new  crack  tip 
(as  shown  in  Figure  62). 


8.2  The  Modified  Willenborg  Model 

When  the  current  plastic  zone  at  the  crack  tip  is  surrounded  by  the 
greatest  prior  plastic  zone,  as  shown  in  Figure  63,  Willenborg,  et  al,  [39] 
assumed  that  the  effective  stress  for  crack  growth  is  reduced  to  a lower 
level  (as  compared  to  the  actually  applied  stress  level)  and  thus  retardation 
in  crack  growth  rate  occurs.  The  load  interaction  zones  ^ and  r^  are 
defined  as  follows: 


and 


r 

P 


F(a)2yT]a 

2 


F(a)2yna0l  /o0L  \2 

o 2 If  ) 

\ ty  / 


(19) 


(20) 


where  F^.  is  the  material  tensile  yield  strength  and  F(a)  is  the  stress 
intensity  geometry  factor,  e.g.,  Equation  4 for  a center  cracked  panel. 

The  parameter  7 reflects  the  state  of  stresses,  i.e.,  7 - 1.0  for  plane 
stresses  and  7=1/3  for  plane  strain;  ^ is  either  equal  to  2.0  (indicates  a 
full  zone)  or  1.0  (indicates  a half  zone,  the  so  called  plastic  zone  radius, 
r^).  Note  Equations  19  and  20  are  the  so-called  Irwin  plastic  zone 
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a CURRENT  CRACK  LENGTH 

r CURRENT  LOAD  INTERACTION  ZONE 

P 

ap  GREATEST  PRIOR  LOAD  INTERACTION  FRONT 
Figure  63.  Load  Interaction  Zone  at  Crack  Tip 
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equations  (i.e.,  same  as  Equation  16  ) . A Rice  equation  can  be  obtained 

just  by  multiplying  a factor  of  1.23374  to  the  Irwin  equation.  As 

previously  discussed  in  Section  7,  the  plastic  zone  sizes  are  most  likely  a 

2 

nonlinear  function  of  (K  /F  ) and  may  also  be  somewhat  affected  by 

max  ty  J 

biaxial  stress  ratios. 


The  stresses  used  in  this  model  are  illustrated  in  Figure  64.  In 
Figure  64 , the  stresses  cr  and  CT  . are  the  spectrum  stresses  in  a stress 
event,  aa^  is  the  applied  stress  which  is  defined  as  an  imaginary  stress 
that  creates  a plastic  zone  whose  boundary  coincides  with  the  greatest  prior 
interface,  i.e. 

F / 2 (a  -a) 

% = F(fy  V Try~~  (21) 

The  stress  CTre(j,  the  reduction  in  applied  stress,  is  defined  as  the 

difference  between  cr  and  the  maximum  stress  in  the  spectrum,  CT  , i.e., 

ap  r max’ 


CT  = a - cr 
red  ap  max 


(22) 


The  effective  stresses  are  then  obtained  as 


(o  ) „ = o - o , 

max  eff  max  red 

(23) 

(CT  . ) rr  = CT  . - cr  , 

mm  eff  min  red 

(24) 

where  both  (o  ) and  (o  . ) rr  are  restricted  to  be  non-negative.  After 
max  eff  min  eff  & 

obtaining  the  effective  stresses,  the  crack  growth  rate  is  determined  from 
the  effective  stress  intensity  range,  effective  stress  intensity  ratio,  and 
the  material  baseline  (constant  amplitude  growth  rate)  data. 
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Figure  64.  Effective  Stresses  in  the  Willenborg  Model 

From  the  above  outline  of  the  model,  the  following  observation  can  be 
made  about  this  model;  (1)  Negative  stress  ratio  is  excluded,  (2)  No 
additional  parameter  is  introduced,  and  (3) 
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If  rr  . > cr  j only  the  reduction  in  stress  ratio  causes  crack  growth 

min  — red  3 

retardation.  Thus  in  selecting  a crack  growth  rate  equation  for  analysis  it 

must  include  a term  to  account  for  the  R effect.  For  <7  . O’  . and 

mm  ^ red 

rr  > (7  ,,  both  Ak  and  R are  reduced.  Finally,  when  <7  < <7  , complete 

crack  arrest  occurred;  this  situation  arises  if  the  overload  ratio  is  2.0  or 
higher.  However,  experimental  data  have  shown  that  an  overload  shutoff  ratio 
of  2.0  is  nonconservative  in  many  cases.  Therefore,  modification  of  the 
model  is  necessary.  Using  a stress  intensity  approach,  Gallagher  & Hughes  [40] 
applied  a correction  factor  to  Willenborg's  reduction  in  applied  stress. 

That  is, 


mod 

rred 


$ cx. 


red 


(27) 


wi  th 


(28) 


This  modification  introduces  two  additional  constants;  (1)  S is  the  overload 
shutoff  ratio,  and  (2)  K^X  is  associated  with  the  fatigue  crack  growth  rate 
threshold  for  R = 0. 


8.3  The  New  Crack  Growth  Retardation  Model 

In  Von  Euw's  dissertation  [34],  he  reported  that  the  crack  growth  rates 
for  the  lower  loads  subsequent  to  a higher  load  were  retarded.  However,  he 
also  observed  that  the  minimum  fatigue  crack  propagation  rate  did  not  occur 
immediately  following  the  application  of  the  overload.  Instead,  the  lowest 
growth  rate  was  reached  after  the  crack  had  extended  over  a distance  within 
the  overload  plastic  zone.  A schematic  representation  of  this  behavior  is 
presented  in  Figure  65.  In  Figure  65,  a ^ and  r^  are  taken  to  be  the  same 
as  those  previously  defined  in  the  Modified  Willenborg  Model  except  that  all 
the  plastic  zones  are  defined  to  be  the  full  plastic  zone  (i.e.,  with  T]  = 2.0); 
a^  is  the  crack  length  after  the  application  of  a high  load  (overload)  and 
r^  is  the  crack  tip  plastic  zone  for  a^  at  a lower  load.  It  was  assumed 
(also  supported  by  experimental  observation)  that  the  crack  would  resume  its 
normal  growth  rate  when  it  travelled  through  the  overload  plastic  zone,  i.e, 
for  its  crack  tip  plastic  zone  to  travel  a distance  of  a*  indicated  in  the 
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top  portion  of  Figure  65.  For  constant  amplitude  loading  with  occasional 
single  overloads,  it  is  convenient  to  assume  that  = r^,  i.e.,  a*  = r^  - r^. 
This  assumption  has  enabled  us  to  construct  the  bottom  part  of  Figure  65. 

In  that  figure,  (da/dN)^  corresponds  to  a^  and  (da/dN^  corresponds  to  a^ 

( = a^  + a*).  The  distance  at  which  the  crack  growth  rate  reaches  minimum 
is  approximately  equal  to  1/8  to  1/4  of  the  total  overload  plastic  zone  as 
determined  by  Von  Euw's  experiments  on  2024-T3  aluminum  alloy. 

To  calculate  the  number  of  delay  cycles  (N  ) for  the  crack  to  travel 
through  a*,  the  following  assumptions  have  been  made. 

1.  Both  (da/dN)^  and  (da/dN^  are  regular  crack  growth  rates,  i.e., 
not  retarded. 

2.  The  crack  growth  rate  behavior  within  the  r^  zone  follows  the 
Paris  Equation  [4l] , and  crack  growth  retardations  are  benefited 
from  "crack  closure." 

3.  (da/dN^.n  occurs  at  x = a*/4 

4.  da/dN  is  linear  between  (da/dN)  and  (da/dN)min,  and  is  also 

linear  between  (da/dN)  . and  (da/dN),,. 

mm  l 

Following  these  assumptions,  we  can  write 

A 

_ dN  1 §Nrni  n . x (29) 

dN  ■ MnJ  la* 

for  0 £ x * £ a , and 

da  _ 1 I, da, 

dN'TRT 

for  £ a < x ^ 


* 

a with  £ = 1/4. 


t (f*£) 
* W2 


(*£)  .A 

aN7  ^dNmi 


(l-|)a 


(30) 
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After  integration,  the  total  number  of  delay  cycles  within  the  overload 
plastic  zone  will  be 


I±L 


''dN,  MN ' . 

1 mm 


(—) 

dNj^_ 


Ahih 


A \A  .A 
W.  VdN'  W. 

mm  Z mm 


(— ) 

J*2 

(— ) 
aN' . 
min 


(31) 


The  resulting  crack  length  versus  cycles  history  is  schematically 
illustrated  in  Figure  66. 

The  crack  opening  concept  is  adopted  to  calculate  the  crack  growth  rates 

inside  the  overload  plastic  zone.  The  effective  Ak  values  in  each  stress 

event  can  be  obtained  following  the  assumptions  set  above  and  the  definitions 

given  in  Figure  67.  There  (AKg)^  = AK^  and  (AK  ^ = AK2«  The  minimum 

effective  AK  would  be  equal  to  the  difference  between  K (K-opening  at 

overload)  and  K at  (a,  + £a*).  Therefore,  this  model  actually  relocates 

the  position  of  (AK  ) . from  a.  to  (a,  + £a“)  and  also  redistributed  the 
v e mm  1 1 

variations  of  AK  from  a small  to  normal  distribution  to  a normal-small- 
e 

normal  distribution. 

Thus,  the  crack  growth  rates  at  various  positions  inside  the  overload 
plastic  zone  are 


(Ak)" 


n 

= C ["(a  -a  . ) J ira1 

max  mm  v 1 


(32) 


C • (AK) 


n 

2 


= C 


(cr  - CT  . ) x/na„ 

max  mm  v 2 


n 


n 


mm 


= c (Ak  ) . 

e mm 


(33) 


(34) 
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Figure  66.  Growth  Rate  Curve  Resulting  from  the  Application  of  a Single  Overload 


and 


The  C and  n are  the  constants  in  the  Paris  Equation  for  the  baseline  material; 
Ci  and  Qgare  Elber's  crack  closure  parameters.  The  crack  growth  rates  between 
these  points  can  be  obtained  by  interpolations. 
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SECTION  IX 


EXPERIMENTS 


A group  of  118  specimens  were  tested  to  investigate  the  effects  of 
biaxial  stress  ratio  on  fatigue  crack  growth  rate  behavior  of  7075-T7351  and 
2024-7351  aluminum  alloys.  The  specimens  were  divided  into  the  following 
groups. 

1.  To  investigate  the  basic  biaxial  ratio  effects  (including  the  applied 
stress  level)  - 57  specimens. 

2.  To  investigate  the  cyclic  stress  ratio  effects  (at  R = 0.7)  ~ 

12  specimens. 

3.  To  investigate  the  behavior  of  cracks  coming  out  from  a hole  - 
10  specimens. 

4.  To  investigate  the  effect  of  biaxial  stress  ratios  on  crack  growth 
retardation  behavior  - 13  specimens. 

5.  To  investigate  the  biaxial  loading  effects  under  out-of-phase 
loading  conditions  - 14  specimens. 

6.  Miscellaneous  tests  including  fatigue  crack  growth  under  sustained 
load  conditions,  angle  cracks  and  fracture  tests  - 12  specimens. 

Tensile  coupons,  center  cracked  panels  (CCT),  and  cruciform  specimens 
were  fabricated  from  10  sheets  of  7075-T7351  and  four  sheets  of  2024-T351 
plate  stocks.  All  the  sheets  of  each  material  were  from  the  same  heat. 

The  size  of  these  commercial  aluminum  plates  was  121.92  cm  (4  feet)  wide  by 
365.76  cm  (12  feet)  long  by  17.7  mm  (0.5  inch)  thick.  The  specimens  were  cut 
from  randomly  selected  areas  in  these  aluminum  plates.  The  testing  conditions 
and  descriptions  of  experimental  procedures  for  each  test  type  are  given  in 
the  following  sections  whenever  is  appropriate.  Discussion  and  correlation 
of  experimental  data  are  also  presented  in  the  same  section  where  testing  of 
a particular  test  type  is  presented. 
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9.1  Material  Characterization  Tests 


Thirty- six  tensile  test  coupons  were  machined  from  all  14  sheets  of 
aluminum  alloys.  Specimen  configuration  was  as  specified  in  ASTM  Standard  E8, 
with  specimen  thickness  equal  to  4.572  mm  (0.18  inch).  A 178,000  N (40  kip) 
MTS  machine  was  used  for  conducting  the  tensile  tests. 

The  average  tensile  yield  strength  (the  0.2  percent  offset  value)  for 
the  7075-T7351  alloy  was  412,022  kPa  (59.8  ksi)  for  both  the  LT  and  the  TL 
directions.  The  average  tensile  yield  strength  for  the  2024-T351  alloy  was 
367,237  kPa  (53.3  ksi)  for  the  LT  direction  and  319,007  kPa  (46.3  ksi)  for 
the  TL  direction.  Engineering  stress-strain  curves  were  also  obtained  from 
each  tensile  test.  One  typical  curve  was  selected  from  each  alloy  and  it 
was  used  for  conducting  the  elastic-plastic  finite  element  analyses. 

These  two  tensile  stress-strain  curves  are  presented  in  Figures  68  and 
69.  A complete  set  of  the  stress-strain  curves,  for  all  36  specimens,  are 
compiled  in  Volume  II  of  this  report.  Tabulated  tensile  properties  for  all 
36  tests  are  also  reported  in  Volume  II  of  this  report. 


9.2  Constant  Amplitude  Biaxial  Crack  Growth  Rate  Tests 

Thirty  7075-T7351  specimens  and  27  2024-T351  specimens  were  tested  under 

various  biaxial  loading  conditions  (-1.5  < o"x/cTy  < 1.75)  at  R = 0.1,  and  at 

various  applied  stress  levels  (0.2  < ff  /F  < 0.6).  The  testing  conditions, 

for  each  test,  are  listed  in  Table  4.  Unless  otherwise  noted,  the  loading 

profiles  for  <r  and  cr  cyclic  stresses  were  such  that  cr  . occurred  at 
r x y } x,  mm 

<y  for  the  negative  biaxial  ratios  whereas  cr  occurred  at  c r 

y,  max  ° x,  max  y,  max 

for  the  positive  biaxial  conditions. 


For  all  tests  involving  compression  loading,  an  antibuckling  guide  was 
used.  The  apparatus  for  preventing  specimen  buckling  consisted  of  two 
square-shaped  aluminum  plates  and  two  circular  steel  plates.  The  steel 
plates  were  inserted  into  the  circular  hole  in  the  center  of  the  aluminum 
plate.  The  test  specimen  was  sandwiched  in  between  the  aluminum  plates. 

The  crack  could  be  seen  from  an  open  slot  19.05  mm  (0.75  inch)  wide  and 
12.7  cm  (5  inches)  long  in  the  center  of  the  circular  plate.  The  circular 
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Specimen  No.  LT-7-2,  E=10,340  ksi 
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Figure  68.  Typical  Tensile  Stress-Strain  Curve,  7075-T7351,  LT,  t = .18  Inch 


Specimen  No.  LT-2-12,  E = 10,390  ksi 


Figure  69.  Typical  Tensile  Stress-Strain  Curve,  2024-T351,  LT,  t=  .18  Inch 


TABLE  4.  BASIC  BIAXIAL  RATIO  EFFECT  TESTS 


Material 

Test 

Case 

,max 

,min 

a 

x ,max 

CTx ,min 

f 

Notes 

7075-T7351 

70 

12.0 

1.2 

-1.8 

-18.0 

8 

7075-T7351 

25 

12.0 

1.2 

-1.2 

-12.0 

10 

G> 

707  5-T7351 

26 

12.0 

1.2 

-1.2 

-12.0 

5 

7075-T7351 

56 

18.0 

1.8 

-1.8 

-18.0 

5 

707  5-T7351 

58 

18.0 

1.8 

-0.9 

-9.0 

5 

7075-T7351 

135 

18.0 

1.8 

-0.9 

-9.0 

2 

707  5-T7351 

64 

30.0 

3.0 

-1.5 

-15.0 

2 

707  5-T7351 

29 

12.0 

1.2 

-0.6 

-6.0 

5 

707  5-T7351 

30 

12.0 

1.2 

-0.6 

-6.0 

10 

707  5-T7351 

27 

12.0 

1.2 

6.0 

0.6 

7 

E> 

7075-T7351 

28 

12.0 

1.2 

6.0 

0.6 

8 

707  5-T7351 

139 

30.0 

3.0 

15.0 

1.5 

2 

7075-T7351 

63 

30.0 

3.0 

15.0 

1.5 

2 

D> 

7075-T7351 

57 

18.0 

1.8 

9.0 

0.9 

5 

7075-T7351 

133 

18.0 

1.8 

18.0 

1.8 

5 

7075-T7351 

55 

18.0 

1.8 

18.0 

1.8 

3 

7075-T7351 

23 

12.0 

1.2 

12.0 

1.2 

10 

G> 

7075-T7351 

24 

12.0 

1.2 

12.0 

1.2 

5 

7075-T7351 

141 

12.0 

1.2 

15.0 

1.5 

5 

7075-T7351 

68 

12.0 

1.2 

18.0 

1.8 

5 

707  5-T7351 

142 

12.0 

1.2 

21.0 

2.1 

5 

7075-T7351 

1 

6.0 

0.6 

0 

0 

5 

E> 

7075-T7351 

2 

10.0 

1.0 

0 

0 

5 

7075-T7351 

123 

4.8 

0.48 

0 

0 

10 

D>  £>  [£> 

7075-T7351 

7 

12.0 

1.2 

0 

0 

5 

7075-T7351 

8 

12.0 

1.2 

0 

0 

10 

s> 

7075-T7351 

125 

4.8 

0.48 

0 

0 

10 

707  5-T7351 

131 

20.0 

2.0 

0 

0 

5 

fl>  E> 
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TABLE  4.  (continued) 


Material 


7075-T7351 

7075-T7351 


2024-T351 

2024-T351 

2024-T351 

2024-T351 

2024-T351 

2024-T351 

2024-T351 

2024-T351 

2024-T351 

2024-T351 

2024-T351 

2024-T351 

2024-T351 

2024-T351 

2024-T351 

2024-T351 

2024-T351 

2024-T351 

2024-T351 


2024-T351 

2024-T351 

2024-T351 

2024-T351 

2024-T351 


y,max 


J . Q U . I - 

y,min  x,max  x,mxn  f 


Notes 


5 E> 

5 Q> 

G> 


2 E> 


8 


b> 


in 


TABLE  4.  (continued) 


Material 

Test 

Case 

°y ,max 

— 

a 

y ,mm 

— 

a 

x ,max 

cr 

x,min 

f 

Notes 

2024-T351 

132 

8.0 

0.8 

0 

0 

5 

[D>  [£> 

2024-T3  51 

14 

10.0 

1.0 

0 

0 

5 

2024-T351 

149 

6.0 

0.6 

0 

0 

10 

C £> 

[T>»  Interferometry  photographs 
[T>  CCT  specimen 

[T>  20Hz  for  da/dN  < 10  ^ inch/cycle 

Precracking  data  for  the  fracture  specimen 
[£>  tl  direction 

[T>  Tested  in  the  biaxial  loading  machine 
[T>  Tested  in  distilled  water 
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plate  could  be  rotated  to  follow  the  crack  growth  direction.  All  the  metal 
plates  in  this  device  were  coated  with  teflon  (approximately  0.06  inch 
thick).  The  outer  plates  were  clamped  to  the  loading  arms  of  the  specimen. 
The  outer  plates  had  set-screw  adjustments  around  the  inner  circle  which  in 
turn  "clamp"  the  circular  plates  to  contact  the  specimen  flat  surface. 
Actually  these  plates  were  allowed  only  to  contact  (i.e.,  the  clamping  load 
was  only  a few  pounds).  From  the  strain-survey  results,  no  noticeable 
effect  of  load  interaction  due  to  the  use  of  such  constraints  was  noticed. 

Both  the  CCT  and  the  cruciform  specimens  were  used  to  develop  crack 
growth  rate  data  for  uniaxial  loading  conditions  (a  = 0).  All  but  two  of 
the  CCT  specimens  were  tested  in  a 356,000  N (80  kip)  MTS  machine.  The  last 
two  CCT  specimens  and  the  cruciform  specimens  were  tested  in  the  biaxial 
loading  frame  to  insure  that  compatible  crack  growth  rate  data  would  be 
developed  from  both  the  CCT  and  the  cruciform  specimens.  The  loading 
characteristics  of  the  newly- built  biaxial  test  unit  were  checked  out 
through  testing  of  the  CCT  specimens  in  both  testing  machines. 

The  size  of  the  CCT  specimens  was  17.78  mm  (7  inches)  wide  by  40.64  cm 
(16  inches)  long,  having  the  central  portion  tapered  down  from  17.7  mm 
(0.5  inch)  to  4.572  mm  (0.18  inch).  The  cruciform  specimen  configuration 
was  discussed  in  the  previous  section. 

Since  in  some  biaxial  loading  cases  the  crack  might  not  grow  perpendic- 
ular to  the  principal  loading  direction,  it  was  desirable  to  measure  the 
crack  in  both  magnitude  and  direction.  A photographic  polar  grid  such  as 
that  shown  in  Figure  5B  was  printed  onto  the  very  finely  polished  cruciform 
specimen  surface.  Tick  marks  were  placed  at  every  15  degrees  around  the 
circumferential  grid  line,  and  the  spacing  between  grid  lines  was  1.27  mm 
(0.05  inch). 
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Both  the  CCT  and  the  cruciform  specimens  were  precracked  from  an  initial 
EDM  slot,  at  the  center  of  the  specimen,  to  the  desired  initial  flaw  size 
(approximately  3.81  mm  in  total  length).  All  specimens  were  precracked  at 
test  maximum  load  level  by  applying  tension- tension  load  cycles  normal  to 
the  EDM  slot  using  the  MTS  machine  or  the  biaxial  loading  frame  whichever  was 
convenient.  Crack  length  measurements  were  made  at  very  small  increments  to 
obtain  an  adequate  understanding  of  the  crack  growth  behavior.  For  the 
cruciform  specimens,  testing  was  terminated  when  the  crack  growth  rate  was 
faster  than  2 . 54  (Um/cyc le  (10  ^ inch/cycle)  or  the  crack  had  reached  the 
border  of  the  flat  area.  In  the  tension-compression  biaxial  ratio  tests, 
testing  was  terminated  when  the  crack  reached  the  end  of  the  opening  slot  in 
the  buckling  prevention  device  (4  to  5 inches  in  total  crack  length). 

The  effect  of  cyclic  frequencies  was  not  the  primary  interest  of  this 
investigation.  However,  due  to  the  nature  of  the  biaxial  loading  tests, 
and  the  operation  characteristics  of  the  biaxial  test  unit,  lower  frequencies 
had  to  be  used  for  testing  at  higher  applied  loads,  whereas  higher 
frequencies  could  be  applied  to  lower  load  test  cases.  Therefore,  as  shown 
in  Table  4 , some  test  cases  consisted  of  several  replications,  and  each  of 
them  was  run  at  a different  cyclic  frequency  to  insure  that  test  results 
would  be  consistent  with  the  range  of  frequencies  being  applied. 

9.2.1  Test  Results 

Stress  intensities  for  cracks  in  the  CCT  specimens  were  computed 
using  Equations  3 and  4.  For  the  cruciform  specimens,  the  test  results 
indicated  that  the  crack  grew  straight  in  all  the  tests  with  ox  < a . 

However,  when  exceeded  a , the  crack  turned  away  from  its  initial  plane 
and  finally  ended  up  propagating  in  a direction  perpendicular  to  ax  (see 
the  example  shown  in  Figure  5B  for  ax  = 1.5  o ). 
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Stress  intensity  values  presented  in  Figure  10  (adjusted  by  the 

actual  0^  in  each  test)  were  used  to  correlate  the  da/dN  data  for  the 

straight  cracks.  A complete  set  of  the  tabulated  a versus  N data  and  da/dN 

versus  K curves  for  each  test  are  compiled  in  Volume  II  of  this  report, 
max 

Composited  da/dN  versus  K curves  for  each  of  the  six  biaxial  ratios  (0, 
r max 

+0.5,  +1.0,  and  -1.5)  are  presented  in  Figures  70  to  79.  Examination  of  all 
the  test  results  revealed  that  all  the  crack  growth  rate  curves  are  almost 
identical,  i.e.,  for  the  same  material  and  cyclic  stress  amplitude,  there 
was  no  effect  on  fatigue  crack  growth  rate  due  to  differing  biaxial  stress 
ratios.  It  is  significant  to  note  that  each  of  these  composited  da/dN 
curves,  (e.g.,  Figures  72  and  76)  are  composed  from  many  sets  of  test  data 
points  and  each  set  of  those  data  had  been  generated  from  different  combina- 
tions of  stress  levels  and  cyclic  frequencies.  Six  test  technicians  were 
involved  at  different  times  in  collecting  the  crack  length  versus  cycles 
records  for  all  the  tests.  Even  so,  it  is  very  evident  that  the  crack 
growth  rate  behavior  for  all  the  tests  remained  consistent.  As  also  shown 
in  Figures  70  to  79,  some  data  points  at  higher  K levels  exhibited  lower 
crack  growth  rate  values.  These  data  points  were  reduced  from  long  cracks 
(a  > 2.75  inches)  in  the  cruciform  specimen.  Extrapolated  K values  were 
used  during  data  reduction.  The  actual  K values  for  these  crack  lengths 
might  have  been  significantly  reduced  due  to  the  reduction  of  stress  at  the 
area  near  the  rim  of  the  specimen. 

The  effect  of  biaxial  loading  conditions  on  crack  growth  can  also  be 
evaluated  from  the  a versus  N curves  plotted  in  Figures  80  and  81.  There 
the  scatter  (on  N)  for  growing  a crack  from  0.5  inch  to  4.0  inches  was 
only  +167,  (327,  total  scatter  band)  for  the  7075-T7351  specimens  and  +177, 

(347,  total  scatter  band)  for  the  2024-T351  specimens.  As  identified  in 
these  figures,  the  scatter  band  for  the  7075-T7351  specimens  at  B = 1.0  was 
already  217,  (as  compared  to  the  327,  total  scatter  band).  For  the  2024-T351 
specimens,  the  347,  total  scatter  band  actually  represented  the  scatter  band 
for  the  -1.0  biaxial  ratio  specimens.  No  definite  trend  in  the  biaxialty 
effects  are  shown  in  these  figures. 
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Figure  72.  Cyclic  Crack  Growth  Rate  Behavior  of  7075-T7351  Aluminum 
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Figure  74.  Cyclic  Crack  Growth  Rate  Behavior  of  7075-T7351  Aluminum 
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Figure  76.  Cyclic  Crack  Growth  Rate  Behavior  of  2024-T351  Aluminum 


Figure  77.  Cyclic  Crack  Growth  Rate  Behavior  of  2024-T351  Aluminum 
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Figure  78.  Cyclic  Crack  Growth  Rate  Behavior  of  2024-T351  Aluminum 
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Figure  79.  Cyclic  Crack  Growth  Rate  Behavior  of  2024-T351  Aluminum 
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Figure  80.  Crack  Growth  Histories  of  the  7075-T7351  Cruciform 
Specimen 
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Figure  81.  Crack  Growth  Histories  of  the  2024-T351  Cruciform  Specimens 


9.2.2  Correlation  of  High  Biaxial  Ratio  Test  Data 

Three  specimens  (Test  Case  Nos.  141,  68  and  142)  were  tested  under 
high  biaxial  ratio  conditions.  The  biaxial  stress  ratios  were  1.25,  1.50 
and  1.75.  The  crack  growth  profiles  for  these  tests  are  shown  in  Figures  82 
to  84.  It  is  seen  in  these  figures  that  the  crack  turned  away  from  its 
initial  plane,  and  adjusted  itself  to  become  normal  to  the  most  dominating 
stress  component,  the  higher  the  biaxial  stress  ratio,  the  sooner  the  crack 
would  complete  the  crack  path  transitions. 

There  were  two  problems  in  the  high  biaxial  stress  ratio  cases.  The 
first  problem  was  how  to  determine  the  K and  the  da/dN  values  and  the 
second  problem  was  how  to  predict  the  shape  of  the  crack  growth  path.  It 
was  apparent  that  the  crack  growth  rate  (under  constant  amplitude  loading 
conditions)  should  be  the  same  for  all  biaxial  ratios  as  long  as  the  K 
values  at  the  crack  tip  were  appropriately  determined.  All  the  three  methods 
discussed  in  Section  6.4  were  evaluated  and  the  results  are  discussed  below. 

The  value  of  da  for  a curved  crack  is  defined  to  be  the  distance 
connecting  two  consecutive  measurements.  Referring  to  Figure  16,  the 
circles  represent  the  polar  photo  grids  printed  on  the  cruciform  specimen. 
Suppose  one  reading  was  taken  at  Point  A and  another  reading  was  subsequently 
taken  at  Point  B;  da  would  be  the  distance  between  Points  A and  B as 
indicated  in  the  figure.  Then  dN  would  be  the  difference  between  the  number 
of  cycles  taken  at  Points  A and  B. 

Stress  intensity  values  for  each  data  point  (raw  data  are  compiled  in 
Volume  II  of  this  report)  were  computed  and  plotted  against  the  associated 
da/dN  values.  The  crack  growth  rate  curves  constructed  by  each  of  the  three 
methods  (for  computing  K)  for  each  of  the  three  tests  are  presented  in 
Figures  85  to  93.  Comparing  all  nine  crack  growth  rate  curves  to  the 
material  baseline  crack  growth  rate  curve  (Figure  94),  it  is  revealed  that 
all  three  methods  generally  correlate  well  with  the  uniaxial  loading  case. 
Table  5 presents  a comparison  of  the  K values  calculated  from  all  three 
methods  for  all  three  high  biaxial  ratio  conditions.  The  comparisons 
were  made  by  superimposing  the  baseline  material  da/dN  versus  K curve  onto 
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Figure  83.  Crack  Growth  Profile  of  7073-T7351 
Cruciform  Specimen  No.  7-28 
Test  Case  No.  68 

Cr  = 18  ksi,  a = 12  ksi,  R = 0.1,  f = 5 Hz 


Figure  84.  Crack  Growth  Profile  of  7075-T7351 
Cruciform  Specimen  No.  7-4 
Test  Case  No.  142 
Ox  = 21  ksi,  Gy  = 12  ksi,  R = 0.1, 
f = 5 Hz 
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Figure  85.  Crack  Growth  Rate  Behavior  of  7075-T7351 
Cruciform  Specimen  No.  7-38 
Test  Case  No.  141 

CT  =15  ksi,  a =12  ksi,  R = 0.1,  f = 5 


da/dN,  inch/cycle 


ksi  / inch 


Figure  86.  Crack  Growth  Rate  Behavior  of  7075-T7351 
Cruciform  Specimen  No.  7-38 
Test  Case  No.  141 

(J  =15  ksi,  CT  =12  ksi,  R = 0.1,  f = 5 

Y 7 \r  7 7 
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Figure  88.  Crack  Growth  Rate  Behavior  of  7075-T7351 
Cruciform  Specimen  No.  7-28 
Test  Case  No.  68 

cx  = 18  ksi,  ay  = 12  ksi,  R = 0.1, 
f = 5 Hz  y 
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Figure  89.  Crack  Growth  Rate  Behavior  of  7075-T7351 
Cruciform  Specimen  No.  7-28 
Test  Case  No.  68 

a = 18  ksi , a = 12  ksi , R = 0.1,  f = 5 
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Figure  90.  Crack  Growth  Rate  Behavior  of  7075-T7351 
Cruciform  Specimen  No.  7-28 
Test  Case  No.  68 

CT  = 18  ksi,  a = 12  ksi,  R = 0.1,  f = 5 
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Figure  91.  Crack  Growth  Rate  Behavior  of  7075-T7351 
Cruciform  Specimen  No.  7-4 
Test  Case  No.  142 

CT  =21  kis,  a =12  ksi,  R = 0.1,  f = 5 


da/dN,  i 


! 

n 
r 


inch 


e Behavior  of  7075-T7351 
en  No.  7-4 
2 


12  ksi , R = 0.1,  f = 5 


1 10  100 
Kmax’  ksi  /I^ 

Figure  94.  Baseline  Crack  Growth  Rate  Behavior 
for  7075-T7351,  LT  Direction  (CCT 
Specimen,  R = 0.1,  f = 5 Hz) 

each  of  these  nine  crack  growth  rate  curves.  The  precision  of  the  calculated 
K values  was  judged  by  comparing  the  K for  the  uniaxial  case  with  a K 
calculated  by  a particular  method,  at  a selected  da/dN  value.  It  is  seen, 
from  Table  5,  that  the  total  sum  of  and  method  seem  to  offer  quite 
reasonable  results.  The  vector  sum  of  and  ^ method  generally  calculates 
a K value  lower  than  expected.  In  contrast,  the  "equivalent  a"  method  some- 
times produces  a K value  slightly  higher  than  expected.  The  perfect  matches 
in  the  upper  and  lower  da/dN  regions  where  the  crack  was  normal  to  and 
respectively,  obtained  from  the  last  method  strongly  suggested  that  stress 
intensity  factor,  crack  growth  rate  and  crack  growth  direction  are  controlled 
by  the  most  dominating  stress,  not  the  principal  stress,  in  a given  stress 
field. 
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TABLE  5.  COMPARISON  BETWEEN  ANALYSIS  METHODS  FOR  A CURVED  CRACK 


METHOD 


Below 

3 x 10-5 

Perfect 

Perfect 

Perfect 

Be  tween 

3 x 10-5 

& -4 

10 

Slightly 

Too  High 

Too  High 

Slightly 

Too  High 

-4 

Above  10 

Too  High 

— 

Perfect 

Perfect 

Below 

3 x 10  ^ Perfect  Perfect 


Be  tween 
3 x 10-5 

&^-4  Too  Low  Perfect 


-4 

Above  10 


Perfect 


Too  Low 


Too  Low 


Below 

3 x 10'5 

Perf ec  t 

Perfect 

Perfect 

Between 

3 x 10"5 

-4 

& 10 

Too  High 

Perfect 

Perfect 

-4 

Above  10 

Too  High 

Perfect 

Perfect 
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A few  words  are  given  here  to  explain  why  the  calculated  K values 
were  too  high  in  the  1.25  biaxial  ratio  case.  As  it  can  be  observed  in 
Figures  82  to  84,  the  cracks  in  the  other  two  higher  biaxial  ratio  tests 
were  relatively  short  and  more  importantly  the  entire  crack  path  of  these 
cracks  were  very  close  to  the  X and  Y planes  of  the  specimen.  Therefore, 
the  magnitude  of  the  stress  components,  c?x  and  a^.,  along  the  crack,  were 
very  much  the  same  as  those  previously  reported  in  Sections  4 and  5.  Since 
the  crack  was  long  and  far  away  from  either  the  X and  Y axis  in  the  1.25 
biaxial  ratio  specimen,  in  the  region  where  the  crack  was  long,  the 
magnitude  of  the  stress  components  (or  the  biaxial  ratio)  might  not  be 
exactly  the  same,  as  could  be  expected. 


After  the  best  suited  procedure  for  calculating  K is  selected,  the 
second  problem  would  be  how  to  set  up  a mathematical  procedure  which  can 
predict  the  entire  crack  propagation  path  of  a crack  starting  from  any 
initial  length  and  angle.  The  computerized  procedure  discussed  in 
Section  8.1  was  applied  to  predict  the  crack  growth  history  for  Test  Case 
No.  142  (B  = 1.75).  In  this  analysis,  the  Collipriest  Equation,  42,  was 
used  to  represent  the  baseline  material  crack  growth  rate  behavior. 

||  = C(Kcr  • AKo)n/2  • EXP [y  • jen(Kcr/AKo)n/s] 

wi  th 

jHAK2/((1  - R)Kcr  • AKq)] 

v = arctanh  

Jta[(l  - R)Kcr/AKo]  (36) 

The  emprirical  constants  Kcr,  AKq,  C,  and  n,  for  the  CCT  specimen  (in  the  LT 
direction)  for  both  aluminum  alloys  were  determined  and  are  presented  in 
Figures  95  and  96.  Numerous  computer  runs  were  made  using  different 
combinations  in  empirical  crack  growth  rate  constants  (i.e.,  different  crack 
growth  rate  curves  were  fitted  between  the  middle  and  the  upper  bounds  of 
the  CCT  da/dN  data).  The  best  two  of  the  computer  runs  are  presented  in 
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Tables  6 and  7.  The  results  are  compared  to  the  actual  test  results  shown 
in  Table  8.  In  these  two  predictions,  one  of  them  exhibited  good  correla- 
tion on  a (too  fast  on  a ) whereas  the  other  one  showed  good  correlation 
x y & 

on  a (too  slow  on  a ) . The  crack  growth  rate  constants  used  in  these  two 
y x 

analyses  were  almost  identical  (only  a small  difference  in  C,  0.016  versus 
0.0155,  while  the  other  three  constants  were  the  same  in  both  analyses). 
Therefore,  it  is  apparent  that  sensitivity  is  very  high  in  this  type  of 
analysis.  However,  it  is  anticipated  that  improvements  can  be  accomplished 
through  obtaining  an  improved  crack  angle  relationship,  or  using  separate 
sets  of  da/dN  constants  for  the  da^  and  the  da^  components,  etc.  An 
important  point  to  be  made  here  is  that  the  present  crack  growth  concept, 
and  approach,  can  be  useful;  a refined  computational  procedure  can  be 
developed  following  the  aforementioned  modifications. 

9 . 3 The  Effect  of  Cyclic  Stress  Ratio 

Twelve  cruciform  specimens  were  tested  at  cyclic  stress  ratio  of  0.7. 
The  testing  conditions  for  these  specimens  are  listed  in  Table  9.  A 
complete  set  of  the  crack  growth  rate  curves  and  the  tabulation  of  a versus 
N data  are  compiled  in  Volume  II  of  this  report.  Typical  crack  growth  rate 
curves  (for  a = 0 and  CT  = + a ) are  presented  in  Figures  97  tolOO. 
Comparing  these  curves,  it  is  seen  that  the  biaxial  stress  effect  is  also 
negligible  at  R = 0.7.  Furthermore,  the  cracks  grew  straight  in  all  the 
tests  (see  Volume  II). 

9 . 4 Crack  Coming  Out  From  a Hole  Tests 

Ten  specimens  of  the  7075-T7351  aluminum  alloy  containing  either  a 

0.25  inch  diameter  hole  or  a 0.75  inch  diameter  hole  were  tested  under  three 

biaxial  stress  ratios  (a  = 0,  a = + 0 , see  Table  10).  The  cyclic  stress 

amplitude  for  all  the  ten  tests  was  0.1.  Two  of  the  ten  tests  were  crack 

initiation  tests.  The  remaining  eight  specimens  contained  very  small  saw 

cuts  at  both  sides  of  the  hole.  Crack  growth  history  for  each  test  was 

recorded.  A complete  set  of  a versus  N data  and  da/dN  versus  K curves 

max 

are  compiled  in  Volume  II  of  this  report.  The  K values  were  calculated 
using  the  superposition  technique  (previously  discussed  in  Section  6.3). 
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TABLE  6.  PREDICTED  CRACK  GROWTH  HISTORY  OF  A 7075-T7351  SPECIMEN 
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Test  Case  No.  142,  B = 1.75,  R = 0.1 
Collipriest' s Crack  Growth  Rate  Constants: 

Kc  = 60.0  ksi  yj  inch  C = 0.016 

AK  = 3.0  ksi  J inch  n = 2.65 


TABLE  7.  PREDICTED  CRACK  GROWTH  HISTORY  OF  A 7075-T7351  SPECIMEN 
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Test  Case  No.  142,  B = 1.75,  R = 0.1 
Collipriest’ s Crack  Growth  Rate  Constants: 

Kcr  = 60.0  ksi  ,/inch  C = 0.0155 

AK  = 3.0  ksi  / inch  n = 2.65 


TABLE  8.  ACTUAL  CRACK  GROWTH  HISTORY  OF  A 7075-T7351  CRUCIFORM  SPECIMEN 
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Test  Case  No.  142.  B = 1.75,  R = 0.1 


TABLE  9.  LOW  AMPLITUDE  TESTS 


Material 

Test 

Case 

CTy }max 

<7 

y ,min 

a 

x ,max 

^x.min 

f 

Notes 

707  5-T7351 

47 

12.0 

8.4 

12.0 

8.4 

5 

707  5-T7351 

48 

12.0 

8.4 

-8.4 

-12.0 

5 

7075-T7351 

49 

12.0 

8.4 

6.0 

4.2 

10 

7075-T7351 

50 

12.0 

8.4 

-4.2 

-8.4 

8 

707  5-T7351 

11 

12.0 
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0 
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G> 
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0 
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7.0 

10.0 

7.0 

10 

2024-T351 

52 

10.0 
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-7.0 

-10.0 

10 

2024-T351 

53 

10.0 

7.0 

5.0 

3.5 

10 

2024-T351 
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10.0 

7.0 

-3.5 
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SPECIMEN  7-51  TEST  CASE  11 


Figure  97.  Cyclic  Crack  Growth  Behavior  of  7075-T7351  Cruciform  Specimen 


SPECIMEN  7-105  TEST  CASE  48 


aXoXo/ipuT'NP/^P 
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Figure  98.  Cyclic  Crack  Growth  Behavior  of  7075-T7351  Cruciform  Specimen 


SPECIMEN  2-19  TEST 


aio/O/ipui  ‘mp/bp 


154 


Figure  100.  Cyclic  Crack  Growth  Behavior  of  2024-T351  Cruciform  Specimen 


TABLE  10 


CRACKS  AT  A CIRCULAR  HOLE  TESTS 


Material 

Test 

Case 

a 

y,max 

O 

y,min 

a 

x,max 

a 

x,min 

f 

Notes 

7075-T7351 

86 

12.0 

1.2 

0 
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5 

7075-T7351 
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0 

0 
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7075-T7351 
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B>  E>  E> 

7075-T7351 

19 

18.0 
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1.2 

12.0 

1.2 

5 

E> 

7075-T7351 

90 

12.0 

1.2 

12.0 

1.2 

5 

7075-T7351 

88 

12.0 

1.2 

-1.2 

-12.0 

5 

•E> 

7075-T7351 

91 

12.0 

1.2 

-1.2 

-12.0 

5 

7075-T7351 

92 

12.0 

1.2 

-1.2 

-12.0 

5 

0.25  inch  di  ameter  hole 
0.75  inch  diameter  hole 
[T>»  CCT  specimen 
[T>  Crack  initiation  test 
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Comparing  the  da/dN  versus  K curves  for  these  cracks  at  a hole  tests  to 
r ° max 

da/dN  versus  K curves  for  specimens  without  a hole,  it  was  observed  that 
max 

the  crack  growth  rate  behavior  was  the  same  in  both  cases.  Again,  the  crack 
in  all  these  specimens  grew  straight  essentially  paralled  to  the  X-axis  of 
the  specimen. 

9.5  Variable  Amplitude  Tests 

Table  11  is  a listing  of  the  testing  conditions  for  13  specimens. 

Single  overload  was  periodically  applied  to  nine  specimens.  The  remaining 
four  specimens  were  subjected  to  spectrum  loadings.  Excepting  Test  Case 
No.  102,  all  the  overloads  had  a magnitude  twice  that  of  the  regular  cyclic 
stress  level  (i.e.,  20  ksi  for  the  2024-T351  specimens  and  24  ksi  for  the 
7075-T7351  specimens).  In  Specimen  No.  LT-2-5  (Test  Case  No.  102),  the 
stress  levels  were  16.7  ksi  for  the  first  two  overloads  and  20  ksi  for  the 
subsequent  overloads.  The  biaxial  ratios  for  these  specimens  were  either 
0,  +0.5  or  -0.5.  The  spectrum  load  tests  were  conducted  at  four  biaxial 
ratios  (0,  -.267,  +.5).  Raw  data  for  all  these  tests  are  compiled  in 
Volume  II  of  this  report.  The  results  are  discussed  in  the  following. 

9.5,1  Single  Overload  Tests 

Crack  growth  history  and  crack  growth  rate  data  for  the  periodically 
overloaded  specimens  are  presented  in  Figures  101  toll.8.  In  part  (b)  of 
Figures  110  to  118,  the  basic  da/dN  behavior  associated  with  a particular  test 
is  revealed  by  excluding  the  retarded  da/dN  data  points.  Comparing  these 
da/dN  curves  with  the  da/dN  curves  for  tests  without  overload  (i.e., 
those  curves  previously  presented  in  Section  9.2),  it  is  seen  that  the 
basic  crack  growth  rate  behavior  in  the  overloaded  specimen  is  the  same  as 
those  in  the  regular  (not  overloaded)  specimens. 
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TABLE  11.  VARIABLE  AMPLITUDE  TESTS 


Material 

Test 

Case 

°y ,max 

or 

y ,min 

(T 

x,max 

c r 

x,min 

f 

Notes 

7075-T7351 

99 

12.0 

1.2 

0 

0 

10 

D>  t£> 

7075-T7351 

103 

12.0 

1.2 

6.0 

0.6 

10 

E>  E> 

7075-T7351 

104 

12.0 

1.2 

6.0 

0.6 

10 

E> 

7075-T7351 

105 

12.0 

1.2 

-0.6 

-6.0 

10 

0> 

7075-T7351 

114 

30.0 

- 

0 

0 

B> 

7075-T7351 

115 

30.0 

- 

- 

-15.0 

E> 

7075-T7351 

116 

30.0 

- 

15.0 

- 

B> 

7075-T7351 

200 

30.0 

- 

- 

-8.0 

2024-T351 

101 

10.0 

1.0 

0 

0 

10 

E>[ l> 

2024-T351 

102 

10.0 

1.0 

0 

0 

10 

IX>D> 

2024-T351 

108 

10.0 

1.0 

5.0 

0.5 

10 

0>[D> 

2024-T351 

109 

10.0 

1.0 

5.0 

0.5 

10 

E> 

2024-T351 

110 

10.0 

1.0 

-0.5 

-5.0 

10 

[D>  Periodic  single  overload,  overload  ratio  = 2.0  (in  both 
X and  Y directions  when  applicable) 

[|^>  Interferometry  photographs 

[5^>  Spectrum  load.  Cyclic  frequencies  vary  (15Hz  for  low 
loads  and  2Hz  for  high  loads) 

[£>•  Tested  in  uniaxial  loading  machine 

[?^>  Periodic  single  overload,  overload  ratio  = 1.67,  2.0 

[!^>  CCT  specimen 
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Modified  Willenborg  Model,  S = 2.3 
^Polynomial  No.  C-15 

E 00 


CYCLES  E 85 

Figure  104.  Crack  Growth  History  of  7075-T7351 
Cruciform  Specimen  No.  7-53 
Test  Case  No.  105 
crx  = -6  ksi,  CTy  = 12  ksi,  R = 0.1, 
f = 10  Hz;  Overload  Ratio  = 2.0 
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Modified  Willenborg  Model,  S = 2.3 
O Polynomial  No.  B-7 
0 Polynomial  No.  B-8 

E 09 


CYCLES  E 06 

Figure  107.  Crack  Growth  History  of  2024-T351 
Cruciform  Specimen  No.  2-38 
Test  Case  No.  108 
ctx  = 5 ksi,  Gy  = 10  ksi,  R = 0.1, 
f = 10  Hz;  Overload  Ratio  = 2.0 
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Modified  Willenborg  Model,  S = 2.3 
APolynomial  No.  B-7 
A Polynomial  No.  B-8 

E 08 


CYCLES  E 06 

Figure  108,  Crack  Growth  History  of  2024-T351 
Cruciform  Specimen  No.  2-34 
Test  Case  No.  109 
ctx  = 5 ksi,  Oy  = 10  ksi,  R = 0.1, 
f = 10  Hz;  Overload  Ratio  = 2.0 
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Modified  Willenborg  Model,  S = 2.3 
□ Polynomial  No.  B-6 

E 00 


Figure  109.  Crack  Growth  History  of  2024-T351 
Cruciform  Specimen  No.  2-41 
Test  Case  No.  110 

a = -5  ksi,  a = 10  ksi , R = 0.1,  f = 10  Hz 
x y 

Overload  Ratio  = 2.0 
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Figure  110.  Crack  Growth  Rate  Behavior  of  7075-T7351,  Cruciform  Specimen  No.  7-114, 
Test  Case  No.  99,  O = 0,  O'  =12  ksi,  R = 0.1,  f = 10  Hz  Overload  Ratio 
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Figure  111.  Crack  Growth  Rate  Behavior  of  7075-T7351  Cruciform  Specimen  No.  7-25 
Test  Case  No.  103;  CTX  = 6 ksi,  cry  = 12  ksi,  R = 0.1,  f = 10  Hz; 
Overload  Ratio  = 2.0 


Crack  Growth  Rate  Behavior  of  7075-T7351  Cruciform  Specimen  No.  7-113, 

Test  Case  No.  104,  <7  = 6 ksi , CT  = 12  ksi , R = 0. 1 , f = 10  Hz , Overload  Ratio 


£ 

o co 

C ti 

•rt  C 

^ O 

04 
•rl 

co  co 

JC  4J 

5 

- Q 

X 

co  ,— c 

S ' 

w <; 


9-[oXo/qout*  np/ep 


o 

0*4 

ii 


171 


Figure  114.  Crack  Growth  Rate  Behavior  of  2024-T351  Cruciform  Specimen  No.  2-24, 

Test  Case  No.  101,  a = 0,  a =10  ksi,  R = 0.1,  f = 10  Hz,  Overload  Ratio 


Figure  116.  Crack  Growth  Rate  Behavior  of  2024-T351,  Cruciform  Specimen  No.  2-38 

Test  Case  No.  108,  CT  = 5 ksi,  CT  = 10  ksi,  R = 0.1,  f = 10  Hz,  Overload  Ratio 


Growth  Behavior  of  2024-T351,  Cruciform  Specimen  No.  2-41 

lase  No.  110,  a = -5  ksi,  C = 10  ksi,  R = 0.1,  f = 10  Hz,  Overload  Ratio 


Two  approaches  were  taken  to  evaluate  whether  the  crack  growth  rate 
behavior  (after  each  overload)  would  have  been  affected  by  the  variations  of 
biaxial  stress  ratios.  A straightforward  approach  was  to  take  a delta  cycle 
count  (between  each  pair  of  overloads)  in  each  specimen  tested.  Comparison 
of  the  delta  cycles  is  presented  in  Table  12.  The  trend  of  the  data  seems 
to  show  that  the  negative  biaxial  stress  field  exhibits  more  retardation 
than  the  positive  biaxial  stress  ratio  of  -0.5  than  a zero  biaxial  ratio 
(uniaxial),  and  even  more  so  when  comparing  the  delta  cycles  between  the 
-0.5  and  +0.5  biaxial  ratios.  This  phenomenon  is  more  evident  when  the 
crack  was  short.  For  a > 1.0  inch  (2  inches  total  crack  length),  the 
biaxial  stress  effect  was  not  clear. 

The  current  state-of-the-art  fatigue  crack  growth  concepts  recognize 
the  importance  of  crack  tip  plastic  zone  size  on  affecting  the  crack  growth 
retardation  behavior.  Crack  growth  analysis  models  always  build  in  a 
mathematical  function  for  computing  the  crack  tip  plastic  zone  size  at  any 
point  along  the  crack  propagation  path.  As  mentioned  earlier,  the  residual 
plastic  zone  sizes,  as  well  as  the  loaded  plastic  zone  sizes,  may  depend  on 
the  combinations  of  and  o . Therefore,  the  general  implication  is  that 
good  correlation  between  a predicted  crack  growth  history  and  a set  of  actual 
test  data  can  be  obtained  by  incorporating  a properly  determined  crack  tip 
plastic  zone  size  relationship  (as  a function  of  stress  level,  crack 
length,  material  mechanical  properties,  state-of-stresses , i.e.,  plane 
stress  or  plane  strain,  and  the  biaxial  stress  ratio,  etc)  into  a crack 
growth  analysis  model.  The  modified  Willenborg  model  and  the  new  retarda- 
tion model  were  used  to  compute  the  crack  growth  histories  for  the  periodic 
single  overload  specimens.  The  modified  Willenborg  model  was  applied  to 
all  nine  tests  but  the  new  model  was  applied  to  the  2024-T351  specimens  only. 
Due  to  the  nonlinear  behavior  of  the  crack  tip  plastic  zone  sizes  (see 

2 

Section  7.3),  estimated  r values,  as  a nonlinear  function  of  (K  /F  ) , 

p max  ty 

were  used  in  the  calculations.  The  estimates  were  based  on  the  backtracked 
residual  plastic  zone  sizes  and  expressed  in  a form  of  polynomials.  Since 
the  backtracked  data  were  built-in  with  the  actual  testing  conditions,  the 
variables,  F(a),  V/2,  and  y in  Equations  19  and  20  vanished.  Thus,  the 
complexities  in  determining  these  geometric  terms  were  eliminated. 
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TABLE  12.  PERIODIC  SINGLE  OVERLOAD  TEST  RESULTS 


ALLOY 

a 

(inch) 

SPECIMEN 

AND  TESTING  CONDITION 

T.C.  99 

T.C.  103  T.C.  104 

T.C.  105 

<CTx  • 0) 

<CTx  ■ V2> 

(CT  = - CT  /2) 
x y 

7075-T7351 

0.25 

\ AN  = 

342  70 

28740 

40530 

0.50 

\ AN  = 

2 7530 

23420  25210 

32660 

1. 00 

► AN  = 

25540 

23810  23650 

1.50 

T.C.  102 

T.C.  101 

T.C.  108  T.C.  109 

T.C.  110 

(CCT) 

(ax  = 0) 

(cr  = a / 2) 
x y 

(c*  - - y2> 

2024-T351 

0.25( 

■ AN  = 

106370 

169790 

0.50' 

AN  = 

93900 

85680  71590 

(Average  = 78635) 

132130 

1.00; 

1 

AN  = 

113590 

(Average 

106950 
= 110270) 

94600  134830 

(Average=114715) 

122060 

1 . 50  ■ 

i 

AN  = 

165300  177810 

(Average  = 171555) 

2.00 
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Since  only  a limited  amount  of  backtracked  plastic  zone  data  was  obtained 
from  the  present  experiments,  it  would  not  be  practical  to  determine  a least 
square  fit  from  each  of  the  figures  (in  Section  7.3).  Therefore,  numerous 
hand-fitted  curves  were  fed  into  the  computer  programs.  The  best  suited 
polynomial (s)  was  found  from  comparison  of  the  predicted  and  the  actual 
crack  growth  histories.  It  should  be  noted  that  in  all  the  following 
predictions,  only  one  r^  function  was  applied  to  both  the  residual  (after 
overload)  and  the  current  (at  the  subsequent  lower  cyclic  loads)  crack  tip 
plastic  zones. 

In  applying  the  modified  Willenborg  model,  the  geometric  parameters 

y and  V in  Equation  21  were  fixed  at  1.0  (for  plane  stress)  and  2.0  (for  full 

plastic  zone),  respectively.  A best  suited  overload  shut  off  ratio,  S,  was 

determined  by  comparing  the  results  obtained  from  numerous  computer  runs. 

Different  combinations  of  S and  r - functions  were  inserted  into  each 

P 

computer  run.  It  was  found  that  an  S value  of  2.3  would  be  best  suited  to 
both  the  7075-T7351  and  the  2024-T351  aluminum  alloys. 

Seventeen  r functions  for  the  7075-T7351  material  were  fed  into  the 
P 

computer  program.  After  evaluation  of  all  the  results,  it  was  determined 
that  Polynominal  No.  C-15  was  best  suited  to  all  five  tests  in  three  biaxial 
ratios  (0,  +0.5).  The  analytically  determined  data  points  are  plotted  in 
Figures  101  to  104  to  compare  with  the  actual  crack  growth  history.  In  each 
test  case,  the  analysis  started  from  the  initial  crack  length  of  that 
particular  test;  each  data  point  plotted  in  the  figure  indicates  the  crack 
length  corresponding  to  a load  Cycle  just  before  the  application  of  an 
overload.  Notice  that  the  first  data  point  always  correlated  best  with  the 
test  result.  No  overload  in  that  period  imp-lies  that  the  Collipriest 
constants  for  the  material  were  properly  determined  and  that  the  baseline 
da/dN  curve  was  accurately  reused.  The  case  of  analytic  data  points 
terminating  in  the  middle  of  a crack  growth  history  curve  implies  that  no 
correlations  were  obtained  beyond  that  point;  either  complete  arrest  or 
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rapid  crack  growth  has  occurred.  It  is  interesting  to  note  that  Irwin's 
plastic  zone  equation  also  exhibited  good  correlation  with  the  actual  test 
results  for  Test  Case  No.  99,  the  uniaxial  case  (see  Figure  101).  The  r^ 
curves  for  the  Irwin  equation  and  Polynomial  No.  C-15  are  presented  in 
Figure  119  for  comparison. 

Another  seventeen  r functions  were  used  to  evaluate  the  plastic  zone 

P 

size  and  biaxial  loading  effects  (using  the  modified  Willenborg  model)  on 
the  2024-T351  material.  In  these  computer  runs,  the  initial  crack  lengths 
were  the  experimental  crack  length  at  the  first  overload  (i.e.,  the  loading 
cycles  prior  to  the  first  overload  were  deleted  from  the  analysis).  After 
examining  all  the  analytical  results,  it  was  noticed  that  the  analyses  were 
extremely  sensitive  to  changes  in  plastic  zone  functions.  Good  correlation 
was  almost  impossible  to  obtain.  Four  polynomials  were  selected  from  all 
the  analyses.  The  analytic  data  points  are  plotted  along  with  the 
experimental  crack  growth  curves  in  Figures  105  to  109.  Polynomial  No.  B-9 
was  considered  to  be  the  best  for  the  uniaxial  loading  condition,  whereas 
the  polynomials  No.  B-7  and  No.  B-8  were  selected  for  the  0.5  biaxial 
loading  condition.  As  shown  in  Figure  109  the  best  selected  polynomial 
(No.  B-6),  for  B = 0.5,  actually  poorly  correlated  with  the  experimental 
data.  Examination  of  the  curves  plotted  in  Figure  120  for  these  plastic 
zone  functions  rather  surprisingly  revealed  that  these  four  functions  are 
almost  identical  especially  between  No.  B-6  and  No.  B-7.  The  reason(s)  for 
such  sensitivities  in  using  the  modified  Willenborg  model  on  the  2024-T351 
material  cannot  be  detected  at  the  present  time.  However,  it  is  certain 
that  considerable  amount  of  modifications  to  this  mathematical  model  will 
be  required. 
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fp,  B-6  = 0.002  + 0.368  <Kmax  /Fty)  - 0.029  + 0.164  (Kmax/F,y)3 

rp,  8-7 . 0.003  ♦ 0.352  (Kmax/F,y)  + 0.004  iKmax/Fty)2  + 0.146  <Kmax/Fty)3 
rp,  B-8  = 0.003  ♦ 0.358  <Kmax/Fty)  + 0.017  (Kmax/F,y)2  + 0.105  IKmax/Fty)3 
rp.  B-9  - 0.002  ♦ 0.356  <Kmax/Fty)  ♦ 0.047  (K^/F,/  + 0.067  (K^/F,/ 


Figure  120.  Crack  Tip  Plastic  Zone  Sizes  in  2024-T351 
Cruciform  Specimens 
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The  2024-T351  specimens  were  reanalyzed  by  using  the  new  retardation 
model.  Again,  various  combinations  of  the  Elber  closure  parameters  (the 
values  for  and  C2)  and  r^  - polynomials  were  fed  into  the  computer 
program.  The  best  correlated  r^  - polynomials  are  presented  in  Figure  121 
and  the  comparisons  between  analytical  and  experimental  data  are  presented 
in  Figures  122  to  126.  In  these  analyses,  the  following  empirical  constants 
were  used. 


Paris  C = 
Paris  n = 

C1 

c2 

T1 

r 


0.01  Micro-inch/cycle/(ksi  J inch) 
2.8265738 
0.5 
0.4 
2.0 
1.0 


The  initial  crack  length  in  each  test  case  was  again  started  from  the  crack 
length  corresponding  to  the  first  overload.  It  can  be  seen  in  Figures  122 
to  126  that  the  correlations  are  much  better  than  those  obtained  from  the 
modified  Willenborg  model.  It  is  also  shown  in  these  figures  that  certain 
biaxial  ratio  goes  with  certain  r^  - polynomial(s)  thereby  indicating  the 
biaxial  effect  on  crack  growth  retardation  behavior. 

In  this  program,  the  last  overload  in  a test  was  often  applied  at 
2 

(K^/F^)  >0.75.  In  the  lower  region,  the  correlations  between  test 

results  and  mathematical  models  were  quite  good;  however,  difficulties 
occurred  at  high  levels.  A literature  survey  was  conducted  of  five 
Air  Force  reports  ^43  to  47]  which  contained  a vast  amount  of  experimental 
data  on  overload  effects.  The  majority  of  the  data  covered  a (K^/F^y) 
range  between  0.1  and  0.25,  and  never  exceeded  0.3.  Therefore,  the 
difficulties  encountered  in  the  high  region  have  never  been  experienced 
in  the  past. 


182 


•c.  A-3  • - 0.004  - 0.529  (Kmax  F,y-  - 1.822  (Kmax  Fty)2  * 5 107  ,Kmax  F,/  -2.665  IKmax  F,y  >4 
'p.  A-7  - -0.004  + 0.465  ;Kmix  Fty.  * 0.017  «max  'F,/ 
rp  A-8  =-0.016  - 0.634  (Kmax  F ty  I + 0.994  ,Kmax  Fly  I2 
rp.R|CE  = r «Kmax  Fty)2  B.O 


Figure  121.  Crack  Tip  Plastic  Zone  Sizes  in  2024-T351 
Cruciform  Specimens 
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Figure  125.  Crack  Growth  History  of  2024-T351 
Cruciform  Specimen  No.  2-34 
Test  Case  No.  109 
CTX  = 5 ksi,  cy  = 10  ksi,  R = 0.1, 
f = 10  Hz:  Overload  Ratio  = 2.0 
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Figure  126.  Crack  Growth  History  of  2024-T351 
Cruciform  Specimen  No.  2-41 
Test  Case  No,  110 

tfx  = -5  ksi,  Gy  = 10  ksi,  R = 0.1, 
..  = 10  Hz;  Overload  Ratio  = 2.0 
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9.5.2  Spectrum  Load  Tests 

Spectrum  load  tests  were  conducted  on  7075-T7351  specimens  at  four 
biaxial  ratios  (0,  -0.267,  +0.5,  see  Table  11).  One  specimen  (B  = -0.267) 

was  first  tested  to  obtain  information  on  the  effects  of  spectrum 
truncations.  During  the  test  duration,  a full  spectrum  (see  Table  13)  and 
a truncated  spectrum  (by  deleting  some  loading  steps  (see  Table  14)  were 
alternatively  applied  and  crack  growth  rate  at  different  test  durations 
were  compared.  It  was  determined  that  loading  steps  containing  stress 
levels  below  one-third  of  the  peak  stress  level  could  be  truncated  without 
causing  significant  distortions  in  crack  growth  rate  behavior.  The 
truncated  spectrum  presented  in  Table  14  was  then  used  for  testing  of  the 
other  three  specimens.  Crack  growth  histories  ("a"  versus  number  of  blocks) 
for  all  four  specimens  are  presented  in  Figure  127.  Crack  growth  rate  per 
block  (da/dB)  versus  K at  the  peak  of  each  block  and  the  crack  growth 
profile  for  each  test  are  presented  in  Figures  128  through  131.  It  is 
evident,  from  the  crack  growth  histories  or  the  da/dB  curves,  that  the 
biaxial  stress  ratios  affected  the  crack  growth  rate  behavior.  Comparing 
the  biaxial  stress  ratios  affected  the  crack  growth  rate  behavior.  Comparing 
the  crack  growth  curves  in  Figure  127  with  the  delta  cycle  trend  reported  in 
Table  12,  it  is  noticed  that  both  single  overload  and  spectrum  load  tests 
exhibited  the  same  biaxial  stress  ratio  trend  on  crack  growth  retardation; 
i.e.,  cracks  grew  faster  at  positive  biaxial  stress  states  but  slower  at 
negative  biaxial  loading  conditions. 

To  analyze  the  spectrum  load  test  results  by  using  the  modified 
Willenborg  model,  it  would  be  logical  to  input  the  previously  selected 
Tp  - polynomial(s) , from  the  single  overload  tests,  into  the  computer 
program.  As  reported  earlier,  plastic  zone  function  No.  C-15  seemed  to  fit 
equally  well  to  all  the  biaxial  ratios  between  +0.5  and  -0.5.  Since  ou  was 
the  same  and  the  loading  steps  were  the  same  in  all  the  spectrum  load 
specimens,  there  was  only  one  analytical  "crack  length  versus  block"  curve 
for  all  four  biaxial  loading  conditions.  Needless  to  say,  it  will  be 
meaningless  to  compare  one  analytic  curve  against  four  experimental  curves 
unless  one  could  prove  that  the  deviations  in  the  test  results  were 
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TABLE  13.  BLOCK  SPECTRUM 


Loading  Step 

P [D> 

max 

p . B> 

mm 

Cycles 

1 

84.92 

5.62 

3 

2 

30.16 

0 

32 

3 

29.52 

7.41 

811 

4 

54.60 

7.41 

103 

5 

28.10 

11.67 

1234 

6 

77.78 

7.41 

6 

7 

44.60 

11.67 

372 

8 

91.75 

5.62 

1 

9 

24.60 

7.41 

778 

10 

68.25 

7.41 

20 

11 

38.89 

11.67 

811 

12 

58 . 2 5 | 

0 

203 

13 

33.33  | 

11.67 

1166 

14 

32.54 

16.03 

264 

15 

35.40 

7.41 

656 

16 

61.90 

16.03 

11 

17 

49.84 

11.67 

279 

18 

60.63 

7.41 

36 

19 

22.06 

11.67 

1834 

20 

20.00 

7.41 

914 

21 

47.78 

0 

2 

22 

79.68 

16.03 

1 

23 

41.43 

7.41 

421 

24 

100.00 

5.62 

1 

25 

12.92 

0 

111 

26 

60.79 

11.67 

97 

27 

23.49 

0 

412 

28 

47.46 

16.03 

61 

29 

54.12 

11.67 

139 

30 

37.46 

0 

8 

31 

74.29 

11.67 

7 

32 

19.21 

0 

1414 

33 

48.10 

7.41 

343 

[T>  Percent  of  peak  stress  in  the  spectrum 
(100  percent  = 30  ksi) 
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TABLE  14.  TRUNCATED  SPECTRUM 


Loading  Step 

p [D> 

max 

p . [£> 

mm 

Cycles 

1 

84.92 

5.62 

3 

2 

54.60 

7.41 

103 

3 

77.78 

7.41 

6 

4 

44.60 

11.67 

372 

5 

91.75 

5.62 

1 

6 

68.25 

7.41 

20 

7 

38.89 

11.67 

811 

8 

58.25 

0 

203 

9 

35.40 

7.41 

656 

10 

61.90 

16.03 

11 

11 

49.84 

11.67 

279 

12 

60.63 

7.41 

36 

13 

47.78 

0 

2 

14 

79.68 

16.03 

1 

15 

41.43 

7.41 

421 

16 

100.00 

5.62 

1 

17 

60.79 

11.67 

97 

18 

47.46 

16.03 

61 

19 

54.12 

11.67 

139 

20 

37.46 

0 

8 

21 

74.29 

11.67 

7 

22 

48.10 

7.41 

343 

E>  Percent  of  peak  stress  in  the  spectrum 
(100  percent  = 30  ksi) 
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7075-T7351  Aluminum,  cr  = 30  ksi 

y max 


Figure  127.  Crack  Growth  Histories  of  Cruciform 
Specimens  under  Biaxial-Spectrum 
Load  Conditions 
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1-INCH  SPACING 


a^oXo/qouf  ‘gp/ep 
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Figure  130.  Spectrum  Crack  Growth  Test  of  7075-T7351,  Cruciform  Specimen  No.  7-63,  Test  Case  No.  200, 
ffy,  max  = 30  ksi,  ^x/^y  = -0.267 


mainly  experimental  scatters,  which  was  apparently  not  the  case  here.  Just 
for  the  sake  of  comparison,  a computer  run  was  performed  using  the  truncated 
spectrum  and  plastic  zone  polynomial  No.  C-15.  The  Irwin  equation  was 
also  used  for  baseline  comparison.  The  analytic  crack  growth  histories  are 
presented,  along  with  the  experimental  data,  in  Figure  127. 

At  this  point,  it  seems  necessary  to  review  the  experimental  data  and 
find  out  why  only  r^  - function  was  obtained  from  the  single  overload 
data  (for  the  7075-T7351  material),  whereas  the  spectrum  load  data  indicated 
the  need  for  four  different  plastic  zone  functions.  The  statistics  on  the 
variations  of  delta  block  counts  between  +0.5  biaxial  ratios  was  extracted 
from  the  raw  data  of  the  spectrum  load  tests  and  are  summarized  in  Table  15. 
Comparing  the  spectrum  load  test  results  (Table  15)  and  the  single  overload 
test  results  for  the  7075-T7351  aluminum  alloy  (Table  12),  it  is  seen  that 
the  variations  in  delta  cycles  between  jK).5  biaxial  ratios  were  41  percent 
and  34  percent,  respectively,  for  the  first  two  pairs  of  overloads  (for  the 
crack  in  these  specimens  to  grow  from  0.5  inch  to  1.0  inch  and  from  1.0  inch 
to  2.0  inches).  However,  the  variations  in  delta  blocks  in  the  spectrum  load 
tests  were  as  high  as  100  percent  in  the  0.5  to  1.0  inch  crack  length  dura- 
tion and  an  average  of  61  percent  in  the  remaining  crack  growth  histories. 
Therefore,  it  is  apparent  that  the  biaxial  ratio  effects  were  magnified  in 
the  spectrum  loading  conditions.  Since  this  experimental  program  only 
included  a few  tests  in  the  variable  amplitude  loading  conditions,  more  tests 
are  required  to  generate  sufficient  information.  Prior  to  that,  development 
of  a useful  biaxial-crack-growth-retardation  analysis  procedure  will  be 
impossible.  However,  it  is  important  to  note  that  the  test  results  generated 
from  this  program  clearly  indicated  that  there  is  a biaxial  loading  effect  on 
crack  growth  retardation  behavior. 
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TABLE  15.  SPECTRUM  LOAD  TEST  RESULTS- 


2a 

(Inch) 

Number  Of  Blocks 

Variations 
In  Ab 
(%) 

B 

= 0.5  B =-0.5 

0.5 

0 

°" 

[•  AB  = 50 

K Ab  = ioo 

100 

1.0 

50; 

100' 

85 

125 

47 

2.0 

135  J 

225; 

► 55 

90 

64 

3.0 

190  j 

3 1 5 J 

► 35 

61 

74 

4.0 

225 

376 

*7075-T735l  cruciform  specimens 
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9.6  Out-Of-Phase  Loading  Tests 

Examination  of  Equations  1 and  2 revealed  that  the  load-stress  relation- 
ship characteristics  in  the  current  cruciform  specimen  design  were  not  the 
same  as  those  in  an  infinite  sheet.  In  the  out-of-phase  cyclic  loading 
conditions,  the  input  loads  for  the  cruciform  specimen  were  adjusted  in 
order  to  make  the  stresses  in  the  cruciform  specimen  at  any  time  during  a 
cycle  compatible  with  the  stresses  in  an  infinite  sheet,  particularity  the 
instantaneous  biaxial  stresses  ratios.  A simple  analytic  procedure  was 
formulated  and  described  as  follows. 

1.  In  an  ideal  specimen,  e.g.,  an  infinite  sheet,  let  o be  the 
reference  stresses;  i.e.,  let  = o°  where  the  superscripts 
i and  o denote  in-phase  and  out-of-phase , respectively.  Therefore, 
the  sinusoidal  waves  for  the  biaxial  stress  components  can  be 
defined  as 


a ° = cr  , mean  + (a  , mean  - cr  , min)  • Sin  GPt 

y y y y 


(37) 


CTX°  = ctx>  mean  + (Ox>  mean  - CTx,  min)  • Sin  (Cbt  +Gi)  (38) 


where  Cdt  is  a time  parameter,  a is  the  phase  angle  and 

1 / , . V 

O , mean  = m (O’  , max  + o , mm) 
y Zy  y 

O , mean  = 4 (o  , max  + o , min) 
x’  2 y’  y 

with  B being  the  biaxial  stress  ratio.  The  instantaneous  biaxial 
stress  ratio  at  any  instant,  Wt,  will  simply  be  the  ratio  of 
Equations  38  and  37,  i.e., 


(39) 

(40) 


y' ideal 


(41) 
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2.  In  a cruciform  specimen,  the  stress  and  load  relationship  is  given 
by 


P 


>12 

1! 

[6.55  - 1.73  Bp] 

(42) 

>JcnI 

II 

[6.57  Bp  - 1.75] 

(43) 

where  the  superscript  c stands  for  cruciform  specimen,  P is  the 

load  applied  on  the  y axis  of  the  cruciform  specimen  and  Bp  is  the 

biaxial  load  ratio  (P  /P  ) in  the  cruciform  specimen.  In  order 

to  make  both  cr  and  cr  and  the  instantaneous  biaxial  ratio  in  the 
x y 

cruciform  specimen  equal  to  those  in  an  ideal  specimen,  set 


and 


therefore , 


) 


ideal 


a° 

y 


ideal 


COt 

6.57B  - 1.75 

P 

cot 

6.55  - 1.73  Bw 
P 


and  thus , 


B^ 

P 


6.55B  ^ + 1.75 
a 

cot 

6.57  + 1.73  Bq 


(44) 


(45) 


(46) 


(47) 


is  the  biaxial  load  ratio  required  for  the  cruciform  specimen  under 
out-of-phase  condition  at  any  instant,  cot. 
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3. 


Using  Equations  42,  43  and  44  the  loads, 
to  a given  BpWT  would  be 


Py  and  Px  corresponding 


P ^ = 12c ° / (6. 55  - 1.73  B Wt) 
y y P 

P <*  = B P 

x p y 


(48) 

(49) 


where  a is  given  in  Equation  37.  Therefore,  the  loads  in  the 

cruciform  specimen,  in  any  given  time  ojt,  corresponding  to  a given 

combination  of  a , a . , B,  and  a,  can  be  determined  by 

y,  max  y,  mm  7 

using  these  equations. 


The  results  of  this  analysis  have  indicated  that  the  adjusted  input 

OJT  OJT . 

loads  (P  or  P^  ) for  the  loading  conditions  combining  R = 0.1  and 

B = -0.5  (or  -1.0),  were  actually  oscillating  between  positive  and  negative 

loads.  This  loading  condition  required  cyclically  loading  the  test  coupon 

from  a negative  load  level  (below  zero)  to  a positive  load  level  (above 

zero)  and  then  dropping  the  load  level  back  down  to  negative.  This  type  of 

zero  load  crossing  operation  would  have  rapidly  worn  out  the  bearings  in 

the  biaxial  test  unit.  Therefore,  a test  matrix  consisting  of  a group  of 

7075-T7351  and  2024-T351  specimens  was  planned.  The  specimens  were 

subjected  to  either  (+)  or  (-)  biaxial  ratios  at  R = 0.7,  or  R = 0.1  with 

only  the  (+)  biaxial  ratios.  Table  16  presents  a listing  of  loads  (P 

7 v y,  max 

P . , P and  P . ) for  B = +0.5  and  +1.0,  a = 180°  and 

y,  mm  x,  max  x,  mm  — — 

R = 0.1  and  0.7.  Here,  R was  defined  to  be  the  ratio  of  C . and  C 

Y j 111  X 11  j 9 uld  A 

and  a was  12  ksi  for  the  7075-T7351  specimens  and  10  ksi  for  the 

2024-T351  specimens.  Note  that  Test  Cases  No.  40  and  117  in  Table  16  were 
not  treated  by  load  adjustments.  The  results  of  these  tests  demonstrated 
the  magnitude  of  the  problem  (the  effect  of  load  adjusted  versus  load  not 
adjusted).  Test  results  arfe  presented  in  Volume  II  of  this  report.  Stress 
intensity  factors  for  these  tests  were  computed  using  the  same  procedure  as 


5 
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TABLE  16.  LOADING  CONDITIONS  FOR  180°  OUT-OF-PHASE  CYCLIC  STRESS  TESTS 
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Interferometry  photographs 

Load  not  adjusted  [5^>  f = 10Hz  [T^>  f = 15Hz 


for  the  in-phase  biaxial  tests.  It  was  noted  that  the  crack  growth  rate 
versus  K curve  for  the  load  adjusted  cases  were  identical  to  those  tested 
in  the  in-phase  loading  conditions.  It  was  also  clear  that  the  crack  growth 
rate  curves  for  the  test  without  load  adjustments  did  not  correlate  with 
the  test  data  previously  obtained  from  the  in-phase  loading  conditions. 

Also  indicated  in  Table  16,  the  biaxial  ratio,  in  each  test,  was  not 
constant  throughout  the  test.  It  actually  varied  as  each  component  of  the 
sinusoidal  waves  oscillate  across  each  other  at  different  phase  angles.  For 
those  tests  conducted  in  the  present  program,  the  biaxial  ratio  could  vary 
by  a factor  of  10  in  each  completed  cycle.  The  crack  was  still  growing 
straight  and  the  crack  growth  rates  were  unchanged. 

9.7  Miscellaneous  Tests 

In  this  section,  experimental  results  for  three  particular  types  of 
testsare  presented.  These  testing  types  were:  sustained  load  tests,  angle 
crack  tests,  and  fracture  tests.  The  test  conditions  are  listed  in  Table  17. 

9.7.1  Sustained  Load  Tests 

In  this  type,  ax  was  held  constant  at  a predetermined  level  (matched 

with  a to  form  a desired  biaxial  ratio).  Only  a was  cyclically 

y , max  y 

loaded  to  propagate  the  crack.  Four  specimens  were  tested  at  plus  or  minus 

0.5  biaxial  ratio  with  R = 0.1.  The  test  results  are  compiled  in  Volume  II 

of  this  report.  Examining  the  da/dN  versus  K curves  for  these  curves  along 

with  the  baseline  da/dN  versus  K curves  presented  in  Section  9.2.1,  it  was 

revealed  that  the  sustained  lateral  stress  component  did  not  cause  any  effect 

on  crack  growth  rate  behavior. 
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TABLE  17.  MISCELLANEOUS  TESTS 


Material 

Test 

Case 

CTy ,max 

°y ,min 

°x,inax 

^Xjmin 

f 

Notes 

7075-T7351 

95 

12.0 

1.2 

6.0 

5 

7075-T7351 

96 

12.0 

1.2 

6.0 

5 

G> 

7075-T7351 

97 

12.0 

1.2 

-6.0 

-6.0 

5 

G> 

7075-T7351 

98 

12.0 

1.2 

-6.0 

5 

G> 

7075-T7351 

79 

12.0 

1.2 

0 

0 

5 

d> 

7075-T7351 

81 

12.0 

1.2 

12.0 

1.2 

5 

707  5-T7351 

83 

12.0 

1.2 

H 

-12.0 

5 

7075-T7351 

150 

. 

- 

■ 

- 

■ 

E> 

7075-T7351 

123 

- 

- 

■ 

- 

■ 

[£>  £> 

7075-T7351 

125 

- 

- 

■ 

- 

| 

[£>  [£> 

2024-T351 

127 

- 

- 

■ 

- 

[t> 

2024-T351 

129 

- 

- 

H 

- 

B 

E>  E> 

B> 

D> 


Sustained  load  test 

Angle  crack  test,  crack  line  oriented  at  45°  with  respect 
to  loading  axes  and  sheet  rolling  direction. 

Fracture  test  (monotonically  loaded  to  failure),  tested  in 
uniaxial  loading  machine. 

CCT  specimen 
TL  direction 
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9.7.2  Angle  Crack  Tests 


Three  cruciform  specimens  consisting  of  a 45°  slot  were  subjected  to 

cyclic  stresses  at  different  biaxial  stress  ratios  (0,  +1.0).  The  specimens 

were  precracked  in  the  biaxial  test  frame  at  a biaxial  ratio  of  1.0,  P = 

max 

25  kips,  and  R = 0.1.  This  method  produced  a sharp  fatigue  crack  from  the 
elox  slot  and  the  crack  was  on  the  same  plane  as  the  elox  slot,  i.e.,  45° 
in  respect  to  the  loading  arms. 

When  the  specimen  was  loaded  to  a -1.0  biaxial  stress  condition,  the 
stress  on  the  45°  plane  was  pure  shear.  Theoretically,  the  crack  should 
have  turned  to  a direction  which  is  70.5  degrees  from  the  original  cracked 
plane  and  then  adjust  itself  to  propagate  in  a direction  normal  to  maximum 
K^.  A sketch  presented  in  Figure  132  shows  that  the  crack  did  turn  away 
from  its  original  cracked  plane.  However,  a pair  of  secondary  cracks  deve- 
loped from  the  sharp  corners  of  the  elox  slot.  The  test  was  terminated  when 
the  secondary  cracks  were  discovered  at  approximately  80,000  cycles  of  test- 
ing time.  Note  that  the  original  crack  grew  approximately  0.02  inch  at  each 
end  and  no  "a  versus  N"  record  was  collected. 

The  second  specimen  (Test  Case  No.  81)  was  tested  under  the  one-to- 
one  biaxial  stress  ratio  condition.  As  discussed  earlier,  the  stresses 
in  all  directions  were  identical  and  the  crack  tip  stress  intensities  were 
pure  Model  1.  Therefore,  the  crack  did  grow  straight  and  stayed  on  its 
original  crack  plane  (see  Figure  133).  The  da/dN  versus  K curve  for  this 
specimen  (see  Figure  134)  was  identical  to  the  baseline  crack  growth  rate 
curves  presented  in  Section  9.2.1. 

For  the  zero  biaxial  ratio  condition,  Test  Case  No.  79,  the  stress 
intensity  at  the  tip  of  a 45°  crack  was  one-half  K ^ and  one-half  (see 
Figure  12).  The  initial  turning  angle  at  on-set  of  crack  propagation 
should  have  been  approximately  55°,  clockwise  from  the  original  crack  plane. 
Again,  the  final  crack  propagation  path  should  have  been  normal  to  the  maxi- 
mum direction  (i.e.,  normal  to  ct  , see  the  actual  crack  growth  profile  pre- 
sented in  Figure  135).  The  K values  for  the  test  data  points  were  calculated 
by  using  Figure  12  and  Equation  15.  Comparing  the  crack  growth  rate  curves 
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Figure  133.  Crack  Growth  Profile  of  7075-T7351  Cruciform  Specimen  No.  7-91 
Test  Case  No.  81,  Angle  Crack,  a = +(T  = 12  ksi,  R = 0.1 
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in  Figures  134  and  136,  the  last  several  data  points  for  this  specimen  fall 
right  on  the  crack  growth  rate  curve  of  the  last  specimen  (Test  Case  No.  81). 
However,  the  first  several  data  points  at  lower  K levels  (i.e.,  at  shorter 
crack  lengths)  did  not  match  the  expected  da/dN  values.  This  was  probably 
due  to  the  fact  that  the  crack  was  turning  from  a 45°  position  to  a 0° 
position,  the  actual  K values  in  the  transition  region  were  never  determined. 


9.7.3  Fracture  Tests 

Crack  growth  resistance  curve  (R-curve)  tests  were  run  for  each 
material  in  both  orientation  (LT,  TL) . CCT  specimens  with  a longer  initial 
crack  length  (a  = 1.0  inches)  were  used.  ASTM  ES61  was  used  as  a guide  to 
the  test  procedure.  (Note:  the  specimen  dimensions  are  not  in  the  ASTM 
recommended  range  for  this  test  type.) 

Specimens  were  instrumented  with  back-to-back  C.O.D.  gauge  calibrated 
to  +17o  accuracy.  Antibuckling  guides  were  used  to  insure  good  alignment. 

The  specimens  were  statically  loaded;  X-Y-Y  recorders  were  used  to  monitor 
the  two  load  versus  compliance  traces.  As  the  compliance  changed,  periodic 
partial  unloading  of  the  specimen  was  required  so  that  the  compliance  slopes 
could  be  computed.  Prior  to  each  unloading,  visual  crack  length  measurements 
on  each  of  the  two  crack  tips  were  recorded.  Thus,  at  each  data  point,  load, 
compliance,  and  optical  crack  length  were  recorded.  Testing  was  continued 
to  specimen  failure. 

Data  analysis  was  performed  using  the  following  procedures. 

1.  The  two  crack  tip  measurements  were  averaged  to  determine  the 
"a -optical ." 

2.  aeff  (effective  crack  length)  was  solved  for  using  Equation  10 
of  ASTM  ES61 . 
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K ^ (effective  stress  intensity)  was  calculated  using  aef£ 
above  and  the  K solution  for  CCT  specimens  noted  earlier. 

4.  The  possibility  of  net  section  yielding  was  noted  if  the  nominal 
stress  exceeded  907.  of  the  average  0.27.  offset  yield  strength  of 
corresponding  tensile  tests.  ("?"  on  Figures  137  and  138) 

A cruciform  specimen  (7075-T7351  alloy)  was  also  loaded  statically 
in  an  84  kip  uniaxial  frame.  The  biaxial  stress  ratio  in  this  specimen 
was  -0.267.  The  specimen  containing  a 0.6  inch  crack  (total  length)  was 
initially  loaded  to  77  kip.  No  crack  growth  was  noticed.  The  specimen 
crack  was  then  cyclically  grown  to  a total  length  of  4.2  inches. 

Monotonic  load  was  reapplied  and  again  reached  the  77  kip  level.  Again,  no 
crack  growth  record  was  obtained,  and  the  test  was  terminated. 
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SECTION  X 


GUIDELINES  FOR  FATIGUE  CRACK  GROWTH  ANALYSIS 


A general  fatigue  crack  growth  prediction  analysis  procedure  consists  of 
the  following  elements: 

1.  Geometry 

2.  Stress  field 

3.  Stress  intensity  factor 

4.  Loads  spectrum 

5.  Mathematical  crack  growth  model 

For  a given  structural  configuration,  if  the  geometric  corrections  to 

the  crack  tip  stress  intensities  are  known,  e.g.,  the  width  correction 

factors  for  the  CCT  specimen,  the  first  step  would  be  to  compute  <j  and  <7  , 

x y 

along  an  assumed  crack  plane,  for  an  uncracked  panel.  Then,  stress  intensity 
factors  will  be  those  discussed  in  Section  6.  For  a biaxial  stress  ratio 
greater  than  unity,  all  the  first  three  types  of  the  cracks  discussed  in 
Section  6,  i.e.,  center  crack,  angle  crack,  and  crack  at  a hole,  will  become 
a curved  crack.  In  this  case,  the  methods  discussed  in  Sections  6.4  and 
8.1  apply. 

In  crack  growth  computations,  existing  mathematical  models  were  found  to 
be  imperfect  in  accounting  for  retardation  in  a biaxial  stress  field.  This 
investigation  was  restricted  in  scope  to  only  evaluating  existing  models, 
although  a limited  amount  of  retardation  model  development  was  carried  out. 
Pending  the  development  of  an  improved  mathematical  model,  the  modified 
Willenborg  model  can  be  used  for  crack  growth  analysis  under  spectrum 
loading.  The  plastic  zone  polynomials  used  in  this  report  (Section  9.5)  are 
not  recommended  at  the  present  time  because  it  is  uncertain  whether  the 
nonlinearity  found  to  exist  is  inherent  in  the  polynomials  or  is  a function 
of  the  specimen  geometry.  For  the  time  being,  it  is  recommended  that  the 
Irwin  and  Rice  plastic  zone  size  equations  be  applied,  respectively,  to  the 
7075-T7351  and  2024-T351  materials. 
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SECTION  XI 


RECOMMENDATIONS  FOR  FUTURE  RESEARCH 


1,  More  periodic  single  overload  tests  should  be  conducted.  The  durations 
for  applying  the  overloads  and  the  magnitude  of  the  overload  should  be 
carefully  planned  so  that  a full  curve  for  r^  (as  a function  of  K)  can 
be  developed.  A large  number  of  tests  should  be  conducted  to  cover  a 
wide  range  of  biaxial  ratios.  Sufficient  number  of  duplications  should 
be  included  in  each  loading  condition  so  that  the  r^  curves  can  be 
developed  by  statistical  means.  The  backtracking  technique  for 
determining  the  plastic  zone  size  is  recommended. 


2.  Additional  overload  tests  on  CCT  specimens  of  various  width  dimensions 
are  needed  to  determine  whether  or  not  the  crack  tip  plastic  zone  size 
nonlinearities  observed  in  this  study  are  attributable  to  specimen 
geometry  effects. 

3.  Negative  R cyclic  crack  growth  data  (both  constant  amplitude  and 
spectrum  containing  tension-compression  cycles)  should  be  acquired. 


4.  An  improved  crack  growth  retardation  analysis  procedure  should  be 
developed . 

5.  More  curved  crack  growth  data  should  be  acquired  at  various  high  biaxial 
ratios.  Spectrum  load  tests,  at  the  same  biaxial  ratio  range,  should 
also  be  conducted. 


6.  Improved  analytic  solution  for  the  crack  turning  angle  relationship  for 
curved  crack(s)  should  be  developed.  A generalized  computer  routine  for 
growing  a curved  crack  is  needed. 

7.  More  fracture  data  (at  various  biaxial  ratios)  should  be  acquired  so 
that  a failure  criterion  for  crack(s)  in  a biaxial  stress  field  can  be 
developed. 
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